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Normal Modes in Overmoded Dielectric- 


Lined Circular Waveguide 


By J. W. CARLIN and P. D’AGOSTINO 
(Manuscript received October 30, 1972) 


The propagation constants (loss and phase) are determined for the nor- 
mal modes in dielectric-lined overmoded (guide radius > wavelength) 
circular guide by numerical solutions of the appropriate characteristic 
equations. It 1s shown that the heat loss of the TM pn modes is quite low 
and decreases with increasing frequency in contrast to the increase in heat 
loss predicted by the perturbation theory for thin linings. The primary 
effect of the low-loss TM, modes on a transmission system using the 
TE: mode ts a reduction in the route bend losses. 


I. INTRODUCTION 


Communications systems are currently under development by the 
Bell System and abroad which utilize as the transmission medium 
dielectric-lined circular waveguide (DLG) excited in the low-loss 
circular electric (TEo1) mode. The presence of a thin dielectric liner, 
bonded to the wall of the waveguide, eliminates the undesirable de- 
generacy which exists between the phase constants of the TEo: and 
TMi modes in unlined metallic guide. It also significantly modifies the 
propagation constant (phase and attenuation) and field distribution 
of many other modes compared to their counterparts in unlined 
waveguide. 
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In order to determine the TEo; transmission characteristics of lined 
waveguide, it is necessary to determine the characteristics of the nor- 
mal modes of the structure. The general problem of determination of 
the modes in a hollow, perfectly conducting metallic cylinder with 
a lossless dielectric lining has been of much interest since the initial 
investigations of Bucholz! and Wachowski and Beam.? 

Unger’ * discusses the basic propagation characteristics of a thinly 
lined DLG suitable for use as a wide-band carrier in the millimeter-wave 
region. He determines the phase velocities of the normal modes exactly, 
but uses a first-order perturbation theory, with lining thickness as the 
perturbation parameter, to obtain the heat loss of the normal modes of 
DLG. It can be shown that the actual heat loss characteristics of the 
normal modes of DLG differ significantly from those predicted by the 
simple perturbation theory for a thin lining, in that modes not of the 
circular electric class may have very low heat loss at high frequencies. 
This was discussed briefly in an earlier paper® which treated the loss 
characteristics of circular symmetric modes (TE,, or TM,,) in DLG. 
The present paper is a continuation of that work with a more complete 
discussion of the characteristics of the normal modes in overmoded 
DLG. Representative values for the purely real normal-mode propaga- 
tion constants based on a numerical solution of the eigenvalue equation 
resulting from the classical boundary value problem formulation for 
a lossless DLG are given. The losses were determined using two differ- 
ent methods. The methods (“induced current” and ‘‘wall impedance’’) 
will be explained in a later section of this paper and we will see that the 
results are consistent. The normal-mode loss behavior in DLG was 
found to have a simple physical explanation in terms of the plane wave 
reflection coefficient of an equivalent infinite impedance plane. When- 
ever possible, the algebraic details are omitted or deferred to a suit- 
able appendix. Some plots of the electric field lines of the normal modes 
and the energy density distribution over the guide cross section are also 
given. 


II. BOUNDARY VALUE PROBLEM FORMULATION 


An idealization of the waveguide structure under investigation is 
shown in Fig. 1. It is a perfectly conducting cylinder of radius 6 to 
which a lossless dielectric of relative permittivity ¢«, and thickness ¢ is 
bonded. The interior region (r < a) is to be filled with an inert gas and 
its electrical properties are those of free-space over the frequency range 
of interest. 
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Fig. 1—The idealized lossless model for dielectric-lined guide. 


The fields of a normal mode (TE>, or TM,,) in the structure may be 
expanded in the usual manner in terms of two scalar functions over the 
two regions of the guide cross section. In the empty region of the guide 
(r < a), they are of the form 


Tr = Neds Out) sin po 

T, = NiJp(Xn1) cos pd 
where NV, is a suitably chosen normalization constant. On applying the 
usual boundary conditions at the free-space dielectric (continuity of 
the tangential electric and magnetic fields) and metallic wall (zero 


tangential electric field) boundaries, we obtain a characteristic equa- 
tion which must be solved for the eigenvalue k,. Symbolically, we have 


r<a (1) 


E(kn, f, 6, €r, p, 1) = 0 (2) 
where k, (the only unknown) is defined as 
x = Xn a, (3) 


p is the circumferential order of the mode, and f is the frequency of 
interest. The fields vary along the z axis as e~*~- for an assumed e%! 
time dependence. The propagation constant h, is related to k, by 


hi = kh? — x4 (4) 


where k is the free-space propagation constant. 
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The development of the characteristic equation is discussed by 
Unger*:* and is outlined in Appendix A for reasons of continuity. The 
notation used in this paper is consistent with that of Unger.‘ It should 
be noted that the normal mode fields corresponding to the solution of 
(2) are not strictly transverse electric (TE) or magnetic (TM) as in 
hollow metallic waveguide. (The circularly symmetric modes are pure 
TE or TM.) However, we carry the usual TE or TM nomenclature 
over with the understanding that this is what the field structure tends 
to in the limit of zero lining thickness. 

For a lossless structure, the eigenvalue k, as given by the solution of 
eq. (2) is either pure real or pure imaginary as is the characteristic eq. 
(2). Because of this, it is always possible to obtain a solution of (2) by 
the well-known “bisection”? method on a digital computer. 

On solving (2) for kn, we can then determine the propagation con- 
stant h, for the TE», or TM,, mode of interest. The differential propa- 
gation constant (AB) between the TEyn(TMp5n) mode and the TE 
mode is an important parameter in the estimation of TEo; loss charac- 
teristics. We define AG as 


AB = hi — hos(rad/m). (5) 


Although the model used is an idealization of the actual lossy struc- 
ture, Af as given by (5) is still an accurate representation of the differ- 
ence in propagation constants. We can estimate the attenuation con- 
stant an(Np/m) for a TEy.(TMon) mode by using a simple physical 
approximation. 

The fields in the lossless structure are approximately the same as in 
the actual lossy guide. Hence they may be used as a very accurate first- 
order approximation to obtain the conduction currents in the lossy 
metallic walls and the displacement currents in the dielectric lining. 
The total heat loss (wall losses and dielectric losses) readily follows. 
We shall refer to this as the “‘induced current’’ method in the following 
discussion. The details are given in Appendix A. We now define the 
differential attenuation constant (Aq), which is also needed to deter- 
mine the TEo; transmission characteristics in DLG, as 


Aa = Qn — Qo}. (6) 


In the preceding section we have indicated one well-known method of 
determining the attenuation and propagation constants of the normal 
modes in DLG. There is another approach, first introduced by Unger,* 
which aids greatly in understanding in an intuitive fashion the normal 
mode characteristics of DLG. This is in essence a quasi-optics ap- 
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Fig. 2—The wall-impedance model for dielectric-lined guide. 


proach, as we replace that portion of the DLG in the region r > a of 
Fig. 1 by an impedance wall at r = a as indicated in Fig. 2. The values 
of the impedances Z, and Z,, which approximate the actual boundary 
conditions relating the tangential fields at the free-space dielectric 
boundary, 








(7) 


Z,==— 


7 =<, 
H, 





are given by a simple plane wave analysis as outlined at the end of 
Section III. We can then obtain a characteristic equation for the equiv- 
alent wall-impedance structure of Fig. 2. Again we defer the details to 
Appendix B. The resultant equation has the form 


E(kn, f, Zz Ze; a) = 0. (8) 


Here Z, and Z, are complex quantities for a lossy metal wall with a 
lossy liner, as is the resulting eigenvalue k,. We use the well-known 
Newton-Raphson method to solve (8) for k, numerically on a digital 
computer. In this case, the propagation constant Y, is complex (the 
fields vary along the guide axis as e~/7 for e##¢ time dependence) and 
given by 

v2, = kh? — x2, (9) 
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where 
ky = Xen 


and k is the free-space propagation constant. We define the differen- 
tial attenuation constant (Ae) and propagation constant (A) as in 
(6) and (7), where it is understood that 


Vn — hn + 1an. 


III. RESULTS 


In this section we give results characteristic of those expected of the 
DLG for use in typical waveguide transmission systems. The results 
are for a polyethylene liner bonded to copper waveguide walls. The 
general behavior exhibited, however, is representative of that to be 
expected of any low-loss dielectric liner on a good conductor. 

Figure 3 is a typical plot of the eigenvalue (k,) for the first four circu- 
lar symmetric modes for the lossless guide of Fig. 1 based on numerical 
solution of the exact characteristic equation. The eigenvalues exhibit 
a cyclic or periodic behavior as the lining thickness increases. The 
eigenvalue for a TE on(TM.») mode in DLG is equal to that of a 
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Fig. 3—Eigenvalues vs lining thickness for the TE,, and TM,» modes, lossless 
guide; wall-impedance model and the exact solution. 
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Fig. 4—Eigenvalues vs lining thickness for the TEi, and TM,, modes, lossless 
guide; wall-impedance model and the exact solution. 


TE. .n—-1(TMo,-1) in unlined guide when the lining is “equivalent” to 
a half wavelength in thickness. The ‘‘equivalent wavelength” in the 
dielectric is given by \/ve, — 1, where ) is the free-space wavelength. 
The interested reader is referred to the discussion preceding (11) for 
details. Similarly, the eigenvalue for a TE..(TM.,) mode in DLG is 
equal to that for a TM,,,.(TE.,»_1) in unlined guide when the lining is 
a quarter-wavelength thick. 

It can be seen that the eigenvalues of the TM,,, modes change rapidly 
with an initial increase in lining thickness for thin linings, while those 
of the TE., mode do not. This is due to a strong electric field at the 
wall for TM,» modes and a vanishing electric field for the TE., modes. 

The TMo: mode eigenvalue behaves in an interesting fashion for 
very thin linings. It decreases rapidly to zero as the lining thickness in- 
creases and then becomes imaginary. This indicates that the TMo: 
mode propagates as a slow wave or as a surface wave closely bound to 
the lining region in DLG. The other circular symmetric modes ex- 
hibit the same surface wave behavior for thicker linings, i.e., a quarter- 
wave lining, half-wave lining, etc. 

The hybrid mode eigenvalues exhibit the same sort of repetitive be- 
havior as seen in Fig. 4. The eigenvalues repeat as the lining thickness 
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increases by a half wavelength. The TM,, eigenvalues change rapidly 
as the lining thickness increases for thin linings and then undergo little 
further change until the lining is a quarter-wave thick. The TE, 
eigenvalue drops to zero and then becomes imaginary for very thin 
linings, indicative of a surface-wave-type mode. The TE,, mode eigen- 
values (n > 1) change very little for thin linings but undergo rapid 
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Fig. 5—Heat loss vs lining thickness for (a) the TE., and TMon modes; induced 
current method and (b) the TM.» modes; induced current method. 
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changes in the quarter-wavelength lining-thickness region. Figures 3 
and 4 also furnish a good comparison of the exact solution (the curves 
shown in Figs. 3 and 4) for the eigenvalues and the ‘‘wall impedance”’ 
solution (theO’s shown in the figures at selected points). The agreement 
is quite good with observed errors of approximately 2 percent or smaller 
for the cases shown. 

In Figs. 5 and 6, the behavior of the total heat loss for some typical 
normal modes of DLG versus lining thickness is shown. We see that 
the TE.» modes are low loss for 0, \/2, \, --- thick linings while the 
TM.n modes exhibit low loss for \/4, 3/4, --- thick linings.’ The 
normal-mode heat loss exhibits the same periodic behavior vs lining 
thickness as the eigenvalues. The TE; loss increases slowly with lining 
thickness for a thin lining; as the lining thickness approaches a quarter 
wavelength, the TE; loss increases rapidly as the mode propagates 
as a surface wave. 

The TEpo: loss for a \/2 lining is less than in unlined guide but still 
greater than that for the TEo1 mode in unlined guide, as it is equivalent 
to the TEo: mode in a slightly smaller unlined guide. 


FREQUENCY = 
100 GHz 


ATTENUATION IN dB/km 


51i-mm DIAMETER 
GUIDE 


o= 5.8 x 10? y/m 
€, = 2.34, 
TAN & = 83 x 1076 





0 500 1000 1500 2000 2500 
LINING THICKNESS IN MICRONS 


Fig. 6—Heat loss vs lining thickness for the TEi, and TMi, modes; induced 
current method. 
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The TMo: heat loss as shown in Fig. 5a increases rapidly for even 
very thin linings as expected for a surface wave mode. The TMon, 
n #1 (Figs. 5a and 5b), heat loss initially increases with lining thick- 
ness, and the initial slope of the loss curve agrees with the first-order 
perturbation theory of Unger,’ but the loss then levels off (+25 um 
thick lining) and falls rapidly to a minimum for a quarter-wave lining 
(650 um). The behavior of the TE,, and TM,, heat loss is similar to 
that of the TM., and TE,, modes as shown in Fig. 6. The initial slope 
of the loss curves for the TE,, and TM, modes versus thickness agrees 
with the perturbation theory predictions also, but from Figs. 5 and 6 
it is obvious that the first-order perturbation theory is valid for pre- 
dicting the normal-mode heat loss for only very thin linings (< 20 um). 
It is also quite clear that many modes (i.e., TEoi1, TEo2, TEy2, TMs, 
TMi) as shown in Figs. 5 and 6 exhibit low heat loss (<5 dB/km for 
a 200-um lining) in DLG. 

Figure 7 is a comparison of the copper wall losses only for TE,, and 
TM.» modes as predicted by the induced current (the curves as shown) 
and the ‘‘wall impedance” method (the O’s shown in the figure at 
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Fig. 7—Copper-loss comparison, wall impedance vs induced current method for 
the TEs, and TM.» modes. 
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Fig. 8—Copper-loss comparison, wall impedance vs induced current method for 
the TEs, and TMi modes. 


selected points). The two methods are in excellent agreement over the 
range of thickness shown. Even in the regions where the losses are 
changing rapidly, the error is only on the order of 5 percent. The 
dielectric losses (not shown) are much smaller than the copper losses, 
but they also agreed to within 5 percent for the two methods. 

Figure 8 is a comparison of the copper wall losses and dielectric 
losses predicted by the induced-current and wall-impedance methods 
for the TE,. and TM,; modes. The two methods differ by approxi- 
mately 3 percent or less for the TMi; mode, but the TEx: difference is 
on the order of 10 to 20 percent for a 150 to 200-um lining. The reason 
for this difference lies in the determination of the boundary conditions 
Z, and Z, in (7). The approximate values of Z, and Z, were obtained 
by use of a plane wave analysis as in Section III. This assumes a mode 
may be represented in the region local to the walls of the guide by 
a superposition of suitably polarized plane waves propagating in the 
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Fig. 9—TE: electric field lines and power density vs lining thickness. 
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Fig. 10O—TM:: electric field lines and power density vs lining thickness. 


DIELECTRIC-LINED WAVEGUIDE 465 


51-mm DIAMETER GUIDE 
a= 58 x 107 u/m 

€,= 2.34, TAN 8 = 83 x 107% 
COPPER WALLS 


AB,RAD/m 





40 50 60 70 80 90 100 110 120 130 140 
FREQUENCY IN GHz 


Fig. 11—TMhi: Ag vs frequency, wall-impedance model. 


the TE,, modes (n ¥ 1), as they have a much lower loss than the 
TM>, modes for thin linings. 

Figures 9 and 10 are plots of the electric field lines and power dis- 
tribution over a quadrant of the guide cross section for the TE. and 
TMi, modes in the free-space region of the guide for various lining 
thicknesses at 100 GHz. 

It can be seen that the lining drastically alters the field lines and 
energy distribution from those in copper waveguide. We also see that 
the low-loss regions of Fig. 6 are consistent with the energy distribu- 
tions shown in Figs. 9 and 10 (as they correspond to the cases in Figs. 
9 and 10 where the energy at the free-space dielectric interface is quite 
low). The TE, energy density also changes much more slowly with 
lining thickness than does that for the TMi: mode as expected from 
the eigenvalue plots of Fig. 4. Note that the TM), field lines and energy 
density change little over a 100- to 600-um range of lining thickness. 
There would be a marked change for a lining thicker than 600 um, as 
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the TM: mode then assumes a surface wave behavior. The energy 
density distribution dependence on lining thickness was found to be 
similar for the circular symmetric TE., and TM,, modes. 

In the preceding discussion, we have considered the normal mode 
characteristics in DLG for unrealistically thick linings. (In Fig. 5, 
the TEo: heat loss is 21 dB/km at 100 GHz for the 600-um thick poly- 
ethylene lining considered.) In the following paragraphs we will 
present some data for more representative lining thickness over a 
typical frequency band (40 to 110 GHz). 

It is well known that a dielectric lining is required to break the de- 
generacy between the TEo; and TM: modes in unlined guide. In Fig. 
11, A@rm,, is plotted for several lining thicknesses versus frequency. 
We see that for a lining thickness of greater than 150 um the increase 
in A@ is relatively small over the frequency range 70 to 110 GHz. We 
also note that A@ is approximately inversely proportional to frequency 
for frequencies greater than 60 GHz for a 200-ym lining. 
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Fig. 12—TEo: copper loss vs frequency; wall-impedance model. 
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Conversely, in Figs. 12 and 13 it can be seen that a lining of 200 um 
increases the TEo: copper loss by 0.12 dB/km at 110 GHz. The dielec- 
tric loss is negligible (~0.01 dB/km) for a 200-ym or thinner lining with 
a tan 6 of 83 X 10-®. On the other hand, the same liner with a tan 6 of 
10-3 would have a significant (0.12 dB/km) dielectric loss at 110 GHz. 
The preceding indicates a polyethylene liner with a thickness of 100 to 
200 um, and a tan 6 < 10~* is suitable for a transmission system using 
51-mm diameter guide. A more precise determination of the optimal 
lining thickness is beyond the scope of this paper, as it depends in a 
complex fashion on the geometrical aberrations of the guide in situ. 

Figures 14 through 16 are A@ plots for several TE and TM modes 
over the 40- to 110-GHz band in 51-mm guide for 0, 100, and 200 um 
polyethylene linings. The TE,, and TE,,(n > 1) mode A@’s are quite 
close to the unlined guide A@’s for a 200-um lining and are essentially 
inversely proportional to frequency. The TM,, mode Ag’s are strongly 
dependent on lining thickness for linings of 100 um or less and essen- 
tially little changed for linings greater than 200 um. The TM,, mode 
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Fig. 13—TEno dielectric losses vs frequency; wall-impedance model. 
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Ap’s shown exhibit the 1/f-type frequency dependence for a 200-u4m 
lining for frequencies > 90 GHz. It should be noted that the TMa 
mode has a very low A@ for a 200-um lining which necessitates the 
avoidance of long mechanical wavelength elliptical distortions over the 
guide cross section. 

The TE,: Af’s are strongly dependent on lining thickness over the 
entire 40- to 110-GHz band, and their behavior is quite different from 
that of the other modes. This is due to the surface wave character of 
these modes. As a consequence of this, the TE,; modes having eigen- 
values >3.83171 (the TEo: eigenvalue in unlined guide) may become 
degenerate with the TEo: mode at selected frequencies within the band 
of interest as shown in Fig. 17. The TEs: mode goes through a degener- 
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Fig. 15—TEi, and TM,, Ag vs frequency; wall-impedance model. 


acy at a frequency <40 GHz, while the TE2; and TE,; modes are never 
degenerate. 

In Figs. 18 through 23, the heat loss of several normal modes is 
plotted for a 0-, 100-, and 200-um polyethylene lining in 51-mm guide 
over the 40- to 110-GHz band. The characteristics may be summarized 
as: 


(«) The TM,, (except for the TMo: mode) heat loss initially in- 
creases with frequency, levels off, and then decreases by an 
order of magnitude for a 200-um lining from 40 to 110 GHz. 
The presence of the lining substantially reduces the heat loss 
at 110 GHz. The losses are as low as 3 dB/km at 110 GHz. 

(zt) The TE,, losses decrease with frequency, but not in the f-? 
fashion as in unlined guide. The presence of a lining leads to 
a frequency-dependent increase in heat loss over the 40- to 
100-GHz band; approximately 5 percent and 30 percent for 
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5mm DIAMETER GUIDE 
€, = 2.34, TAN 8'=83 x 1076 


AB, RAD/m 
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FREQUENCY IN GHz 


Fig. 16—TEon and TMan A vs frequency ; wall-impedance method. 


the TEo: mode at 40 and 110 GHz, respectively, for a 200-um 
lining. 

(it) The TE»: and TMo; heat loss is quite high (> 100 dB/km) for 
even thin linings, as these are surface-wave-type modes. 

(iv) The TE gn (n #1) heat loss is low (<10 dB/km), and the 
presence of the lining actually reduces the heat loss over much 
of the 40- to 110-GHz range. 

(v) The heat loss is substantially different from that predicted by 
the simple perturbation* theory for thin linings (for linings 
>10 wm) except in the case of the TE., modes. 


The normal mode characteristics, and in particular the low loss of the 
TM,, modes, discussed in the previous paragraphs (Figs. 3 through 
23) can be easily understood on recognizing that the lower-order modes 
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5i-mm DIAMETER GUIDE 
€. = 2.34, TANS= 83 x 1076 


r 


175-~m THICK LINER 


AB, RAD/m 


—20 


—40 
40 50 60 70 80 90 100 
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Fig. 17—TEn: AB vs frequency; exact solution. 


have a plane-wave-type field structure. The electric field of the plane 
waves may be polarized perpendicular (TE... modes) to or parallel 
(TM., modes) to the plane of incidence as shown in Fig. 24. A TE, 
or TM,, mode (for p #0) is necessarily a superposition of both 
polarizations. 

In examining the interaction of these plane waves with the walls of 
the guide, it is helpful to consider the transmission line equivalent 
structures of Fig. 24. The transmission line parameters (characteristic 
impedance Z, and propagation constant k, in the radial direction) are 
dependent on the polarization. In the empty region of the waveguide 
they are given by: 


Perpendicular Polarization Parallel Polarization 


TEon Modes TM on Modes 
k, = k,/a Khe = ha/t 
ka ka (10) 
Z.=n— Ze=n— 


kn ka 
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= 58 x 107 u/m ~~ _=_ a 
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Fig. 18—TM.,, heat loss vs frequency; wall-impedance model. 
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a= 5.8 x 107u/m 
€, = 2.34, TAN 6 =83 x 1076 
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Fig. 19—TE,,» heat loss vs frequency; wall-impedance model. 
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Fig. 20—TEii and TM, heat loss vs frequency; wall-impedance model. 
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Fig. 21—TEi2 and TE; heat loss vs frequency; wall-impedance model. 
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Fig. 22—TE2: and TMon heat loss vs frequency; wall-impedance model. 
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Fig. 23—TEx: and TE,; heat loss vs frequency; wall-impedance model. 
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Fig. 24—Equivalent transmission line structures for the TE.» and TMon modes. 


where 7 is the impedance of free space, k, the eigenvalue for the mode 
of interest, a the guide radius, and k the free-space propagation con- 
stant. In the dielectric region we have for the near grazing incidence 
(6; = 90°) case: 


Parallel Polarization 
TMon Modes 


Perpendicular Polarization 
TE., Modes 





k,&kvVe,—1 k, & kvVe, — 1 ait) 
Ze n/Ver — 1 Z, & wve, — l/er 


The single-layer case may be generalized to the multilayer case by 
simply adding the equivalent transmission line sections. The copper 
walls of the structure have a skin-effect surface impedance on the order 
of 0.05 + 7 0.05 to 0.1 + 70.1 ohms at 40 and 110 GHz, respectively, 
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for both polarizations. On transforming these impedances in the nor- 
mal way through a 200-um polyethylene liner with equivalent trans- 
mission line parameters as given by (11), we obtain the wall impedances 
Z, and Zs which are seen at the dielectric free-space interface. These 
are on the order of 0.05 + 7 50to0.1 + 7 200 ohms. On the other hand, 
the characteristic impedances as given in (10) for the lower-order 
modes of most interest in the empty region of the guide are on the 
order of 10 to 80 ohms for parallel polarization (TM. modes), while 
the perpendicular polarization (TE.,) characteristic impedances are 
on the order of 1000 to 10,000 ohms. 

The reflection coefficients of the dielectric-clad copper wall for the 
equivalent transmission line structures in Fig. 24 are a function of 
Z./Z- OY Z3/Zc. The propagation constant h, or attenuation constant 
a for the mode of interest is related to the phase or amplitude of the 
appropriate reflection coefficient, respectively. In unlined copper guide, 
both of these ratios (Z./Z. and Z4/Z.) are <1. For the TM, modes in 
DLG, |Z.| is on the order of or much greater than Z,, and thus these 
modes are drastically changed from their copper guide equivalents. 
This is shown by the large changes in A for these modes and the drastic 
reduction in the heat loss. On the other hand, |Z 4| is much less than 
Z. for the TE.» modes, and hence the Ag and heat loss of these modes 
is little changed from that in copper guide. The TE,, modes (n = 1, 
p * 0) have a TE,,-like field structure in the region of the wall in 
copper guide and thus are little affected by the lining. The TM,, 
and TE,; modes have a TM.,,-like field structure in the wall region in 
copper guide, and they are drastically altered by the presence of a lining 
as we have observed. 


IV. SYSTEM IMPLICATIONS 


In the preceding sections, we have seen that the TMi and TE: 
modes have very low heat losses ($3 dB/km) at 110 GHz. This low 
heat loss leads to a reduction in route bend losses as can be seen in the 
following. Unger’ shows that the TE: mode undergoes added losses in 
a route bend due to two different effects. The first effect, called mode 
conversion loss, arises from the fact that a portion of the TEo: energy 
is converted into other modes (predominantly TMi and TE, in 
DLG) in the bend and is not reconverted back to the TEo; mode but 
instead emerges from the bend as a spurious mode. 

The second loss mechanism, defined as the added heat loss, is associ- 
ated with the energy initially converted from the TEo: mode to the 
TE. or TMi: mode and then reconverted back to the TEo: mode as 
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the bend is traversed. Since this energy travels in the higher-loss TE.2 
and TM, modes over a portion of the bend, it suffers additional at- 
tenuation over and above the TE: heat loss in DLG. For a given mode, 
the mode conversion loss and added heat loss in a route bend of radius 
R and total bend length L, with a taper length ¢ (the radius of curvature 
is assumed to have a linear taper of 0 to R from 0 to ¢ and RF to 0 for 
L — ¢to L), is given by 











L Mitivion a2 << 
RB—AH (aade eat i0ss) = — Aan 
R? Apz 
(12) 
; 4 (2 
Lrs—mc (mode conversion loss) = . 
(Rf)? ABs 


Here C is the normalized coupling coefficient! due to curvature of the 
guide axis, and the other parameters are as previously defined. In 
Figs. 25 and 26, the total added TE,, loss, due to the route bend, as 


ADDED HEAT 


TEy2 
ADDED HEAT 


5i-mm DIAMETER GUIDE 
200-um LINER 

€, = 2.34 

LOSS TANGENT = 83 x 1076 
5-90° BENDS/km 

61.0m (200 FEET) BEND RADIUS 


TAPER LENGTH 10m 
(53% OF LENGTH IS IN 
BENDS) 


ROUTE BEND LOSS IN dB/km 





40 60 80 100 120 140 
FREQUENCY IN GHz 


Fig. 25—Theoretical route bend losses. 
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IN BENDS) 


TE, 
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ROUTE BEND LOSS IN dB/km 


™ 4, 
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Fig. 26—Theoretical route bend losses. 


well as the individual components comprising it is plotted for some 
idealized routes comprised of five 90-degree bends with a constant 
radius of curvature (61 m) and linear tapers of 1 or 10 m in length in 
DLG over a length of 1 km. 

From Fig. 26, it is obvious that the mode conversion losses are far 
more significant then the added heat loss in route bends with a short 
taper (or zero) length. On the other hand, for a reasonable taper length 
(10 m), the mode conversion losses are significantly lower than the 
added heat loss arising from coupling to the TMj; and TE: modes. 
The linear taper is probably a reasonable approximation to that ex- 
pected in practice as the radius of curvature will change in a smooth 
fashion as the waveguide progresses from a straight run into a curve. 

The added heat loss is substantially reduced by the presence of the 
lining due to the great reduction in the TM; heat loss. For example, 
at 110 GHz the TM, mode yields an added heat loss for the route 
bends of Fig. 25 of approximately 0.005 dB/km. If the TMi; mode 


DIELECTRIC-LINED WAVEGUIDE 479 


heat loss had followed the perturbation’ theory prediction, the added 
heat loss for the same route bend would be ~0.250 dB/km. 


Vv. CONCLUSION 


We have discussed the basic propagation characteristics of the nor- 
mal modes of dielectric-lined overmoded circular waveguide. Some of 
the results obtained differ in surprising fashion from those predicted 
by the perturbation theory of Unger.? The most significant result is 
the low heat loss of the TM, modes at the upper end of the 40- to 
110-GHz band discussed. This effect greatly reduces the added heat 
loss for the TE; mode in route bends over this part of the frequency 
spectrum. 

A simple physical explanation based on an equivalent transmission 
line problem derived from an analogous plane-wave problem was ad- 
vanced to support the results contained here. An experimental investi- 
gation of the decreasing TMi; heat loss* vs frequency has been com- 
pleted, and the data confirm the trends predicted herein. 


VI. ACKNOWLEDGMENT 


The authors are indebted to W. DeLang for his perspicacity and 
thoroughness in carrying out the numerical computations. 


APPENDIX A 
Exact Formulation 


The waveguide structure under consideration is shown in Fig. 1 
of the main text of this paper. The notation used conforms with that of 
Unger.* Much of the following is similar to that contained in Unger?! 
but is repeated here for reasons of continuity and completeness. 

The fields in the guide can be derived from two suitable scalar wave 
functions composed of suitable products of trigonometric and Bessel 
functions. 


Tn = Nad p(Xn1) sin pd 


0O<r<a (13) 
T1, = Nond p(XaP) cos ph 
x2 Loon (Xar 
Th = Nna—J (kn) EO po 
Xe Z pa (KS) 
a<r<b (14) 
x2 Bon (X41) 
Tr = Nn—Jp(kn) ———— os po 
Xe Z pn (ka) 
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where 
ZorlXhr) = HE (Xr) — cHY(xX4r) 
ZorlXar) = HY (xar) — HS? (X4r). 


Zn and Z,, are the derivatives of the above expressions with respect 
to the entire argument. The constants c and c’ will be defined later. The 
T functions satisfy the wave equation 


1fa oT 0 /1 oT 
wea t[2 (2) 4202) - or as 
r Lor or Od \r dd 


where x” is given by x2 or X° for0 <r < aora <r < brespectively, 
and V? is the transverse Laplacian operator. 

The fields in the guide may be derived from the scalar wave functions 
as in Unger‘ 











OT n or; 
E, = > Vn |—"+ a, | 
n or rod 
He=>V [— q oon 
we bred. BF 
OT'n h2, aT, 
n rdg ke Or 
H SI — 24g h3, | (16) 
a "Lor " he rod : 
x2 
Hz = jwee >, Vedu ss Ti, 
‘ (17) 
; x? 
E, = joo Xu sae fe 


where yo and eo are the permeability and permittivity of free-space, 
respectively, ¢ is the relative permittivity (1 for 0 <r<a, e, for 
a <r<_b) over the guide cross section, k = wvuoece is the intrinsic 
plane-wave propagation constant in the various regions of the guide. 
Here, ¢, k, and X have constant but different values over the guide 
cross section as may the separation constant d,. A mode for the struc- 
ture consists of one term in the series given in (16) and (17). For a given 
modal field distribution, the individual V, and J,’s in (16) and (17) 
represent a single forward or backward traveling wave, and they are 
related to each other by the modal impedance Z, 
Va ha 
Z,=—: (18) 


Tx WE 
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Imposition of the appropriate boundary conditions at r = a, andr = b 
leads to a determination of dn, c, c’, Xn, and Xj; on solution of a trans- 
cendental equation. From Unger* we have 


HY (pks) 
° (oka) 
| HY (pks) 
”  HD"(pks) 


(19) 


(20) 


where the ’ associated with the Bessel function denotes differentiation 
with respect to the entire argument. We also define p, kn, and ky, 


p = b/a (21) 
kn = Xn (22) 
ke = X£a. (23) 


We find d, is constant over the guide cross section and is given by 


ae ee a (24) 
dn 7 pk?a*(e, — 1) PaRAnE cae 
where ; 
ade Jp (Kn) 
oe ed (Ba) 
Lipa ns 
Kong (k%) = Se (25) 
KS Zon (KS) 
Zon ks 
¥ on (ki) = Fon (hr) 
ks Bon (ke) 


k, and k; are related by 
he (w €ou0 = h2)a? (26) 
ke? = (wereouo — h2)a?. 


The individual modes may be determined by solution of the eigenvalue 
equation: 


Ly pn(kn) — Fon(ki)]-Lypn(kn) — &r Ypn(ki)] = per — 1)? 


2 [L274 
n 


k4 ke (27) 





for quasi TE,,, TM,, modes, p #0. For circular electric waves 
(p = 0), the eigenvalue equation is 


[Yon(Kn) aetl ¥ on(kS) | = 0. (28) 
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For circular magnetic waves (TMon) we have 
CY on(Kn) — €r Y ales). = 0. (28a) 


The constant NV, is given by imposing an ortho-normalization cri- 
terion on the transverse fields given in (16). This condition is 





n=m 
0 nwAm. 


Valm 





On evaluation of this integral, we find that N, is exactly given by 


T d?, h2, p? 
3 NA kh So (hn) {(2 - x (2 er + kh ypn + 200] 
a Ble) (04) 2 
k3, e,k?/ ksh 
Kn\? dh hi Z3n(pks : 
a (ea wey gerd 
ks] ke? Z3n(k$) pr kz 
p? 
shi tpminest) 


Z3,(pks) p? ) | 
r{ p2 ———_ — } 1 - — ke? V2, OY <. = 1 29 
i (+ fas i pa t RE YB. + OY, (29) 


for TE,, or TMyn modes, p ¥ 0 and by 

















aN? k2 J3(k 





d?, h?, 
2 {a + k2, Yrn + 2Y on) 


+ (2) [oe eet pers dare] = 
ke pP Z3n(kS) n on on 
WN, ki J3(kn) {a + kRyzn + yon) 
ii ; (=) [e Zin(pkn) = (1 ae kee Y2 + QV || = ] 
r ies p Zanlk) n on on 


for circular electric (TEom) and magnetic (TMo,) modes, respectively. 
We shall use a subscript » when referring to TE,, or TM,, modes and 
a subsequent m for TE on, modes. 

Solution of the appropriate eigenvalue equation (27) or (28) leads 
to a value for the propagation constant (hm or hn) which is pure real 
as we have assumed a lossless dielectric-lined guide in the above analy- 
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sis. For an actual guide, h, and h» have a small imaginary part which 
we call the attenuation constant (@om OF @pn) for the mode in question. 
The loss is due to currents flowing in the metal walls of the guide and 
the lossy dielectric. 

The heat loss or attenuation constant (a) due to currents flowing in 
the metal walls of the structure is given by 


_1Pi 


"9 Pr 
where 


1 1 1 211 
pr= > fiseas=— ff (H.-H; +H,-H4)| Rzbdode, 
2 2Jo Jo rab 
i.e., power loss/unit length, 
Pras [8.x Hlds => [Bex H.dS = Vln 


i.e., the total average power flowing through the guide cross section 
for a given mode, and R, is the surface resistance of the metal walls. 





H (b) af adn N.Xx2J (i ) Z pn( pk) cos pd 
z — n&n OF a Ss 
co oie eC aes 
x4, Z'pn( pk’ h2, Z pn(pké, 
xe Digit KS) w? wo €9 €r bb Bpn(ki) 


noting I, = Vi/Kn where Ky = hn/weo. We find 
_ WAR. 


copper 2 





N, n(= =) 3k) 


Zonlokn) mn) dnh2p Z pn(pk§) 
Zon( IS) ~ @otpig ened Z onl KS) 


s = Yi Bon(ok) 
Who Z pn(k§) 
We find for TE), modes using a similar argument that 
_ Uhn dm|\? | Zom(pk)é 
oom * RslNal* (= =) rate (4 ) st el 
copper We ke, Wp9 : 


Zom(k)&, 
and for TMo, 
| Bynlokis) | 


Th,aR Wer €g\?- 
Beene = iv.l?(= "sain ( xs 
copper WE hn ZnlkS) 
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The dielectric heat loss may be found in a similar manner as we now 
have 


] 1 oH -ck* 1 fi 2 ce 
Pr=- [ wear = =| aV=- | | i oH - E*rdrd¢édz 
2 2 o 2Jo Jo a 


where a is the conductivity of the dielectric. Using field approximations 
in the lining as given in Unger we have 


—V,A | (xs S)+ (Ce a) cos (+ — “) ks | sin po 
VA (7) ~ (dn Xs 0) | sin (+ 2 “) ks | eeene 


xX, r : 
E, = jl,AS —sin (+ _ *) ks | sin pd 
a 


WEQ 


A=WN is k 
‘ (=) Tes) 
_ 1 a 


~ gin dks’ ~ Gos ake 
By substituting for the fields and performing the integration the 
power loss is obtained. Note that o = we, eo X [lining loss tangent ] 
where e,; is the relative dielectric constant of the lining. 
For the TE,, or TM,, modes 


Ith, ka\4 
Cnn ssslectric a Ler uoes tangent) ||N.|? (=) J2(kn) 


(6 + 2)6k? ‘ =) | 
ee ees 1 ae 2 a; 2 
x | : | s*(1 +2) + crlas| 


+ 4[ks sin 26k% + 3(1 — cos 26k%) ] 





E, 


By 


where 





xe2 
x | s* (1 _ ~~) C?|d,| | + SC[Re (dn) ]p sin 26k¢ 


nr 





4(S? + C2|d,|2) + 4 sin 26k4(—S? + o*|ds1)]}- 
We find for TE,,, modes that 
= TIh»[e-(loss tangent) ]|Nn|? (= ) 2(km) 


X [(6 + 2)6ke? — kg, sin 26k%, — $(1 — cos 26k%,) ] 


ai ie 





a 
°™ dielectric 
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and for TMo, modes 


a 
°” dielectric 


Ilhn lea\4 
Shae [e-(loss tangent) ]|NVn|2 (=) 6(Kn) 


e2 
x 8? {c 4 2) aks (1 “ =) + [hg sin 2ak¢ 





x92 
+ 3(1 — cos 26k) ] (1 — -) 


i) 


APPENDIX B 
Wall Impedance Formulation 


The determination of the wall impedance (Z, and Zs) are discussed 
in Section III of the text. Here, we will outline the derivation of the 
eigenvalue equation and the subsequent determination of the propaga- 
tion and attenuation constants. The boundary conditions at r = a are 


Z¢ —_ vey ds or 
Lz = = B,/A 5| yaa: 


As before, we can express the fields in the simplified wall impedance 
structure of Fig. 2 as a superposition of two scalar functions. 
Tn = NaJ(Xar) sin po 
Th = Nad p(XnP) Cos pd. 


(30) 


(31) 


Here p is an integer, and we must solve a suitable eigenvalue equation 
to determine X,, as before 


Xn = Ky (382) 
and 7, the longitudinal propagation constant is given by 
X2 = kh? — 2. (33) 


The tangential fields at the wall (r = a) are given by 




















OT n oT n 
Ey = Vn | —d, | a 
rod or r= 
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E, = Jwuo ie ke fie ae 
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On applying the boundary conditions (80) we obtain the desired 
eigenvalue equation 
Xn 
nF alin) | Tolle) + fea Zy~ Soll) | 
j0€o Lg — —— i = 0 
p hi ; Xn 
— FFs) + Tolln) | Tolls) + joes Zs = Soll) | 
where J,(kn) is defined 


d 
(Kn) = 7) 





(35) 


zt=kn 


The solutions (k,) of this equation were determined numerically on 
a digital computer using a Newton-Raphson iterative scheme. The 
wall impedances Z, and Z, are complex for the most general type of 
wall impedance waveguide as is the eigenvalue k,. On determination 
of kn, we obtain the complex propagation constant Y, from (383). The 
attenuation constant a, (propagation constant h,) is then simply the 
imaginary (real) part of Yn. 
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Experimental Verification of Low-Loss TM 
Modes in Dielectric-Lined Waveguide 


By J. W. CARLIN and A. MAIONE 
(Manuscript received October 30, 1972) 


Experimental measurements have been carried out to verify the recent 
theoretical prediction that, in dtelectric-lined circular waveguide, modes 
other than those of the TE on family can exhibit a loss characteristic which 
decreases with increasing frequency. The TM, waveguide was launched 
in 51-mm diameter circular waveguide lined with a 200-micron layer of 
polyethylene, and its loss characteristics were measured. Losses of 23 
dB/km at 80 GHz and 3.6 dB/km at 110 GHz were measured, thus con- 
jirming the decreasing loss with increasing frequency trend predicted by 
theory. 


I. INTRODUCTION 


Recent analytical studies!~* of the heat loss characteristics of the 
normal modes in oversized dielectric-lined circular waveguide indicate 
that modes other than those of the circular electric type have low 
loss (<10 dB/km). Specifically, in 51-mm diameter waveguide lined 
with a 200-um layer of polyethylene, the TM,; mode was predicted to 
have a theoretical attenuation coefficient which decreases with in- 
creasing frequency and achieves values less than 10 dB/km over the 
80- to 110-GHz frequency range. This paper reports some recent ex- 
perimental results which confirm the decreasing heat loss versus fre- 
quency trend predicted for the TM; mode in lined waveguide. 

The differential phase constant (Af) for the TEo: and TM; modes, 


Ag =h, — ho, (1) 


was also measured and found to be in good agreement with the theoreti- 
cal predictions. Here, hn and ho, are the phase constant (propagation 
constant in lossless guide) for the spurious mode and the TE»: mode, 
respectively. 

487 


488 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1973 


II. MEASUREMENT METHOD 


The loss and AS measurements were carried out with a CW re- 
flectometer test set as shown in Fig. 1. Before proceeding to a dis- 
cussion of the test set, it is helpful to first summarize some recently- 
arrived-at conclusions‘ concerning the normal-mode scattering pro- 
perties at the junctions of three common types of circular waveguide. 

The normal modes of helix or dielectric-lined guide, as defined herein, 
are characterized in terms of the equivalent wall impedance of the 
guide. They are defined as TE,, or TM, modes, although they are 
in fact hybrid modes. The so-called TE,, or TM, n modes tend to the 
TEpn or TMpn modes of metallic waveguide as the wall impedance (or 
equivalently the lining thickness for dielectric-lined guide) tends to 
zero. Since the modes in dielectric-lined or helix guide may differ from 
the modes in copper guide, we will add the appropriate superscript 
(Cu-copper guide, DL-dielectric lined guide, HX-helix) to the usual 
TE,, or TM,» notation. 

For the helix and dielectric-lined guides used in the present study, 
the scattering properties of interest over the 80- to 110-GHz frequency 
range may be summarized as follows: 


(i) The TES? mode scatters £80 percent of its energy into the 
TEI mode. Almost all of the remaining energy scatters into 
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Fig. 1—CW reflectometer test set. 
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the TM?* modes which have much higher losses then the 
TEH* model. 

(ii) The TES? mode scatters >90 percent of its energy into the 
TMP? and TER” modes. 

(iii) The TEX* mode scatters >95 percent of its energy into the 
TMV? mode. 


The loss of the waveguide sample is measured using the circuitry 
in the upper portion of Fig. 1. A dominant TED mode is converted 
to a TES mode via a transducer and passed through a smooth copper 
taper to a 51-mm output diameter. At this point, the copper wave- 
guide is connected to a helix waveguide followed by a dielectric-lined 
guide (as shown in Fig. 1) or simply to a dielectric-lined guide. The 
loss of the sample per unit length can then be determined from the 
slope of the reflected signal level versus piston position. A lossy re- 
flecting piston (<—6 dB reflection coefficient) is used to minimize 
the effects of trapped resonances in the test set. The slide screw-tuner 
as shown in Fig. 1 is required to tune out the residual mismatch in the 
system. 

The Cu-HX-DL tandem connection results in the excitation of a 
relatively pure high-loss TE{i* mode which, in turn, scatters almost 
all of its energy into the low-loss TM” mode as indicated by the 
reflectometer trace in Fig. 2. The trace in both the HX and DL sec- 
tions is smooth and free of large oscillations which would be present if 
other modes were strongly excited. The slope of the trace in Fig. 2 
may then be used to determine the loss of the TEX* or TM?" mode. 

On the other hand, the Cu-DL arrangement results in the excita- 
tion of both the TM?" and TER" modes, which leads to the oscillatory 
trace of Fig. 3. The period of the oscillation in Fig. 3 is one-half of the 
beat wavelength (\z) between the TE” and the TMQ” modes, where 
Az is defined as 


\sp = | 2a/ (ho edt = hrm?) te (2) 


h is the phase constant for the modes. For a 51-mm guide with a 
200-um-thick polyethylene lining, the analytic prediction for the 
TMB! — TEP? beat wavelength at 68 GHz is 0.38 m as compared to 
the experimental value of 0.36 m for the trace in Fig. 3. Note that the 
relative magnitude of the oscillation is decaying rapidly as the piston 
travels down the DL section (the TE)” mode has very high loss) and 
thus we could use the slope of the average value of the reflected signal 
in Fig. 3 (the dashed curve) to obtain the TMP loss. 

The test set in Fig. 1 may be modified in a simple fashion to obtain 
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Fig. 2—Returned signal vs piston position for the helix-dielectric guide system, 
TMDL loss measurement. 


the differential propagation constant (Ag) defined in (1) for the TEQ” 
and TMP" modes. On introduction of a small deliberate tilt at the 
Cu-HX interface (as shown in Fig. 1), the TE{? mode excites the 
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Fig. 3—Returned signal vs piston position for the copper-dielectric system, TMP" 
loss measurement. 
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Fig. 4—Returned signal vs piston position for the helix-dielectric guide system 
with a deliberate tilt, TMP" — TED" differential phase Ap. 


TE and TE modes and subsequently the TM?” and TER” 
modes. The beat wavelength (Az) between the TEo” and TM?” modes 
is then twice the period of the oscillations of the moving piston trace 
shown in Fig. 4. The differential propagation constant A@ is now re- 
lated to the beat wavelength Az 


AB = 2a/Xz. (3) 


III. RESULTS 


By making the moving piston runs described above at different fre- 
quencies, it was possible to infer the loss and A@ vs frequency char- 
acteristics of different waveguide samples. Figures 5 and 6 are the loss 
and phase characteristics obtained for the TE? mode in a 20-m-long 
sample of copper waveguide. Here, the 51-mm diameter copper output 
taper can be connected directly to the 20-m sample of 51-mm diameter 
guide. This arrangement excites both the desired TE{? mode and a 
small amount of TM? (—13 dB down) due to TES’ — TM{? mode 
conversion in the taper. It was found that the undesired TM mode 
could be lowered to approximately 15 dB under the TE{? level by 
placing a 5-m length of helix waveguide between the taper and the 
20-m sample of copper guide, and thus this configuration was used. 

The TE{, measured loss appears to be 15 percent larger on the 
average than the theoretical prediction for the TE loss in perfectly 
straight and smooth copper guide. This increase is similar to that ob- 
served in dominant-mode rectangular guide® and to the measured 
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Fig. 5—TE§p loss vs frequency. 


TEo1° copper loss in 51-mm copper waveguide and has been attributed 
to surface imperfections. The individual data points in Fig. 5 are 
accurate to within +10 percent, which accounts for the +4-dB spread 
about the average. The accuracy of the above CW measurement 
technique could be improved if a better taper were available. 
AB(TEY? — TES?) also agrees well with the theoretical value as shown 
in Fig. 6. The experimental points differ by <3 percent from the 
theoretical curve. In summary, the method used yields experimental 
data in good agreement with the theoretical predictions for the at- 
tenuation and differential phase characteristics of the TES? mode. 
Figures 7 and 8 are comparisons of the experimental results for the 
TM?” loss and differential phase constant, A8(TMPY — TED), in 
dielectric-lined guide with the theoretical predictions obtained using 
the methods outlined in Carlin and D’Agostino.? The copper walls 
were assumed to have an effective 15-percent increase in loss due to the 
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Fig. 6—-TE: — TE differential phase Ag vs frequency. 


effects of surface imperfections. The actual waveguide tested had an 
adhesive layer of Hydropol* of thickness 12.5 um to which a 200-um 
layer of polyethylene was bonded. The electrical properties (dielectric 
constant, ¢,, and loss tangent, tan 6) of the Hydropol and polyethylene 
layers were measured’ at 70 GHz and found to be 


Polyethylene Hydropol 
Er 2.28 2.5 
tan 6 10-3 1.4 X 107% 


These values were used in obtaining the theoretical heat loss. A sample 
30 m in length was used in the measurement of the TM?" loss and 
differential propagation constant. 


* Partially hydrogenated polybutediene manufactured by Phillips Petroleum Co. 
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Fig. 7—TMbP? loss vs frequency. 


The qualitative behavior of the TM}! loss, decreasing with increas- 
ing frequency, is in good agreement with the theoretical prediction. 
The measurement accuracy is approximately + 10 percent at the 
lower frequencies (80 to 90 GHz) and +20 percent at the higher fre- 
quencies (100 to 110 GHz). The quantitative agreement is also good 
from 100 to 110 GHz. The discrepancy at the lower frequencies (80 
to 100 GHz) is probably due to the fact that the lining thickness is not 
a uniform 12.5 ym or 200 um for the dielectric layers as assumed in 
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Fig. 8—TMDL — TED? differential phase Ag vs frequency. 


the theoretical models. This has two effects: 


(2) 


(72) 


The total heat loss as shown in Fig. 7 as well as the losses in the 
various layers and the copper walls is strongly dependent on 
the thickness of the dielectric layers. A liner with an “average” 
thickness of 200 um would have an “‘average”’ heat loss higher 
than that expected for a smooth 200-ym liner, as the loss in- 
creased more rapidly for a decrease in lining thickness than vice 
versa, as seen in Fig. 7. 

The field structure of the TM2” mode fluctuates more rapidly 
and interacts with fluctuations in lining thickness to a greater 
degree as the frequency decreases. This is intuitively reasonable 
from the plots of the TM{i" mode energy density over the guide 
cross section as given in Carlin and D’Agostino.’ There the 
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energy density at the air-dielectric interface in dielectric-lined 
guide decreases as frequency is increased for a fixed lining 
thickness, or alternatively decreases as lining thickness is in- 
creased for a fixed frequency. This suggests that mode con- 
version losses for the TMQ” mode arising from lining thickness 
fluctuations increase with decreasing frequency in dielectric- 
lined guide. A more quantitative investigation of this effect is 
beyond the scope of this paper. 


The measured differential phase constant (AQ) is in excellent agreement 
with the theoretical prediction as shown in Fig. 8. This is further con- 
firmation of the fact that the launching arrangement in Fig. 1 does 
indeed excite the TM" mode in dielectric-lined guide. 


II. CONCLUSION 


In summary, it has been confirmed experimentally that the TM)" 
dielectric-lined waveguide mode is a mode which has an attenuation 
constant which decreases with increasing frequency, thus confirming 
qualitatively the behavior predicted by theory.'* The measured 
values (~ 3.75 dB/km) are approximately 30 percent greater than the 
theoretical predictions from 100 to 110 GHz. At 80 GHz, the measured 
value (23 dB/km) is approximately 2.5 times larger than the theoretical 
prediction. The large discrepancy at the lower frequencies is probably 
due to TMRY — TER” mode conversion losses arising from fluctuations 
in the lining thickness. 
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Prefilters, Sampling, and Transmission 
Rates for Intraframe Codecs for 
Picturephone® Service 


By R. T. BOBILIN 
(Manuscript received September 21, 1972) 


This study of prefilters and sampling rates analyzes the edge busyness 
phenomenon which ts a very significant degradation introduced by intra- 
frame video coders. For any given transmission rate, the prefilter and 
sampling rate can be varied to yield a tradeoff between the edge busyness 
and the rise time associated with a video transition. A narrower prefilter 
results in less edge busyness but slower rise times and vice versa. Similarly, 
varying the sampling rate implies an inverse variation in the number of bits 
per sample allowed in the quantizer, which leads to a tradeoff between alias- 
ing and quantizing noise, both of which are components of the edge busy- 
ness phenomenon. This paper shows in detail how the prefilter, sampler, 
and quantizer affect both edge busyness and rise time, and optimizes the 
intraframe coder design as a function of transmission rate from an analysis 
of this two-factor tradeoff. 


I. INTRODUCTION 


The goal of this paper is to study the effects of filtering and sampling 
in the digital processing of signals for Picturephone® service. In doing 
this we are ultimately interested in choosing the optimum prefilters 
and sampling rates (and therefore the number of bits per sample used 
in the quantizer) for intraframe DPCM codecs. To simplify the analy- 
sis it is assumed that the edge busyness phenomenon is the most sig- 
nificant degradation introduced by these codecs. It will be shown how 
this edge busyness is affected by each part of the codec and how it can 
be controlled by trading it off against picture resolution. The main 
theme of this study is this tradeoff between the detrimental effects of 
both edge busyness and degraded rise times. Both factors are associ- 
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ated with video edges and hence only edges, modeled as step functions, 
will be considered. 

Section II is a brief description of the analog Picturephone system, 
which outlines relevant short-haul system specifications. Section III 
defines edge busyness and shows why it is such a significant video im- 
pairment; it also describes how edge busyness is affected by each part 
of the codec and how it can be controlled by trading it against picture 
resolution. This edge busyness—resolution tradeoff is further investi- 
gated in Section IV which uses a series of subjective pair comparison 
tests to determine the optimum tradeoff as a function of system quality. 
Section V uses the results of Sections III and IV to show that optimum 
coder composition, in terms of prefiltering, sampling rate, and quan- 
tizer structure, is a function of the digital transmission rate. This study 
is concluded in Section VI with a subjective testing program which 
substantiates the validity of the edge busyness—rise time evaluation of 
video coders and rates a number of intraframe coders on a five-com- 
ment impairment scale. 


II. ANALOG SYSTEM DESCRIPTION 


For the purposes of this study, the important parts of the analog 
Picturephone system can be modeled as illustrated in Fig. 1. Here both 
the camera system and all analog links connecting the transmitter to 
the receiver are nominally flat to 1 MHz. If all the frequency shaping 
in the receiver (excluding de-emphasis) is combined into one equivalent 
rolloff filter, the resulting filter has a crispened Gaussian response that 
rolls off to —20 dB at 1 MHz. Hence the total optics-to-optics step re- 
sponse of the analog Picturephone system is dominated by that of the 
receiver station set.* This total optics-to-optics rolloff and overshoot 
characteristic is given by :t 


[R(f) |? = {1 + K(f/T)?] exp [—0.5/T)? I}? (1) 


where f is in MHz and: 
K = 0.5292 


T = 0.35592. 


From this the step response can be shown to be (assuming linear phase 
with a slope of to): 


1 Erf (B eis B 2 
s(t) = 3 BB) aD eS ) (2) 


* See Ref. 1, pp. 291-292. 
See Ref. 2, eq. (8), but add a set of missing parentheses. 
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Fig. 1—Analog system model. 


where 
B= v2xrT(t — to) X 10°. 


This step response is plotted in Fig. 2. The response is symmetric about 
the half-amplitude point with a 4-percent single undershoot and over- 
shoot and has a 10- to 90-percent rise time (7'z) of 0.69 us. 


III. DIGITAL SYSTEM DESCRIPTION 


The essential parts of the intraframe DPCM digital system model 
are illustrated in Fig. 3. In addition to the analog system already de- 
scribed, there is a prefilter, and there are sampling circuits and a quan- 
tizer-decoder combination. The prefilter limits line interference and 
shapes the analog video signal to be encoded. The sample-and-hold and 
sampling clock circuits turn the analog video into a discrete time for- 
mat. The quantizer is used to classify the sample-and-hold outputs into 
a finite number of values for transmission over a given digital channel. 
The decoder is the inverse of the DPCM quantizing algorithm and 
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Fig. 2—Analog system step response. 
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Fig. 3—Digital system model. 


puts out the quantized version of the coder’s sample-and-hold wave- 
form to the receiving analog system. 

In order to see what the codec does to a video transition (such as the 
edge of a face or a shirt or any black-white boundary), let us start with 
the simplest part of the codec and study its effect on the video signal. 
Then the remaining parts of the codec can be added one at a time and 
analyzed. All step responses will be measured at the output of the 
analog receiver, as would be seen by the user. 


3.1 Sampling Circuits 


Consider the simplest video codec consisting of only a sample-and- 
hold circuit driven by a sampling clock. Hence the prefilter, quantizer, 
and decoder circuits of Fig. 3 are deleted. If we model the ouput of the 
Picturephone transmitter as a perfect step function (the validity of this 
model will be discussed in the next section), the input to the codec 
sample-and-hold circuit is illustrated in Fig. 4a. If the sampling instants 
of the coder clock are given by sampling Phase 1 in Fig. 4a, then the 
outputs of the sample-and-hold circuit and the receiver rolloff filter are 
those designated Phase 1 in Figs. 4b and 4c, respectively. If the phasing 
of the sampling clock with respect to the step function was changed by 
a quarter of the sampling period, the sampling instants could be given 
by sampling Phase 2 of Fig. 4a. The resulting outputs of the sample- 
and-hold circuit and the receiver rolloff filter would be given by the 
Phase 2 curves in Figs. 4b and 4c, respectively. The outputs from this 
second phasing of the sampling clock would be the same as those of the 
first phasing delayed by 7/4. If the sampling clock were delayed as 
given by sampling Phase 3 and 4, the corresponding outputs would be 
given by the Phase 3 and 4 waveforms in Figs. 4b and 4c, respectively. 
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In this way we can see that the step response of this digital system 
(always observed at the output of the receiver rolloff filter) will always 
have the same shape as that of the analog system. The time of occur- 
rence of the step response will depend on the relative phasing of the 
sampling clock with respect to the camera step function. This indeter- 
minacy of the edge location due to the sampling phase dependence is 
one specific example of aliasing and will henceforth be called sampling- 
induced busyness (SIB). Quantitatively, this sampling-induced busy- 
ness will be measured as the maximum indeterminacy of the edge at 
the 50-percent amplitude point. This measure is illustrated in Fig. 4c 
and in this case is equal to 7’, one sampling interval. 

Sampling-induced busyness is detrimental because it causes the 
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Fig. 4—(a) Step input from camera. (b) Decoder output. (c) Rolloff filter output. 
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breakup of edges and is easily seen in a typical Picturephone display. To 
see how this breakup occurs, turn to Fig. 5. This represents an enlarge- 
ment of a small section of a Picturephone display. Only five (L1 to L5) 
of the 267 lines in a frame are illustrated, and only four sampling in- 
stants (S1 to S4) are shown. The main diagonal represents the edge of 
a black-white transition. In the original scene all the area to the left of 
the edge is black and all area to the right of the edge is white. On line 
L1, the black-white transition is detected by the sampler at sampling 
instant S4. On line L2, the transition is also detected at sampling in- 
stant 84. It is detected at 83 on line L3 and at S2 on both lines L4 and 
L5. Thus the transition is detected at the times indicated by circles on 
Fig. 5. After these sample-and-hold waveforms are passed through the 
receiver rolloff filter, the edge of the black-white transition would ap- 
pear to the observer to be the crooked dashed curve of Fig. 5. 

For a codec whose sampling clock is locked to the raster, the 
above curvature is motionless if the edge is still. If the codec is not 
locked to the raster or if the edge is moving from frame to frame, the 
above curvature will crawl through the picture. A straight line will be 
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Fig. 5—Sampling-induced busyness. 


INTRAFRAME CODECS 503 


turned into a moving crooked edge. This effect is extremely visible 
and quite annoying for an unfiltered sampled-and-held picture. 

The only effect of introducing the sampling circuits into the analog 
system is this edge indeterminance we have called sampling-induced 
busyness. The rise time of the system has not been affected by the 
sampling circuits. 


3.2 Prefilter 


The addition of the prefilter changes the input to the sample-and- 
hold circuit from a perfect step function to a gentler transition with 
a nonzero 10- to 90-percent rise time as might be given in Fig. 6a. The 





2-MHz SAMPLING RATE 





o—> 
4—_e 
— 
—_ 
— 
— 


2T 3T 4T 5T 6T 


2T 3T 4T 5T 6T 


t t t f t t 


2T 3T 4T 5T 6T 7T 8T 
TIME, t —> 


Fig. 6—Effect of prefilter on SIB. (a) Prefilter output. (b) Sample-and-hold output. 
(c) Picturephone® receiver output. 
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Fig. 7—Sampling-induced busyness—rise time tradeoff. 


more prefiltering, the gentler the transition. The effect that this pre- 
filter has on the decoder step response is illustrated in Fig. 6b for four 
particular sampling clock phasings (in particular those phasings given 
byt = 0, 7/4, T/2, and 3/74 on the time scale of Fig. 6a). As opposed 
to the no prefilter case where the step functions are simply delayed in 
time (see Fig. 4b), the addition of the prefilter results in an overlapping 
of the decoder step responses. After these decoder outputs are passed 
through the receiver rolloff filter,* the sharp corners will be smoothed 
and the width of the step response envelope will be reduced from the no 
prefilter case, as can be seen by comparing Fig. 6c with Fig. 4c. Thus 
sampling-induced busyness can be reduced at the expense of degraded 
rise times. The exact tradeoff between sampling-induced busyness 
(SIB) and average rise time (7'z) is given in Fig. 7 for a Gaussian- 

* See the Appendix for an outline of the computer program used to derive system 


step responses for any combination of transition step size, prefilter step response, 
sampling rate, and quantizer. 
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shaped prefilter. The subjective weighting between these two factors is 
dependent upon picture content and will be discussed in Section IV 
after the quantizer effects are added. 

The curve of Fig. 7 corresponds to a fixed sampling rate. Higher 
sampling rates will shift the curve closer to the analog case, and lower 
sampling rates will shift the curve further away from the analog case. 
This effect is illustrated in Fig. 8 for no prefiltering and for a Gaussian 
prefilter 6 dB down at 1 MHz. The higher the sampling rate, the 
closer the knee of the tradeoff curve approaches the analog condition 
of no busyness and a 0.69-us rise time. Other families of prefilter shapes 
(non-Gaussian) produce slightly different SIB-7'2 curves, but all ex- 
hibit the same tradeoff, reducing one only at the expense of increasing 
the other. 

At this point it is necessary to discuss the validity of modeling the 
output of the Picturephone camera system as an ideal step function. 
The camera frequency characteristic is nominally flat to 1 MHz; 
therefore, as explained in Section II, the optics-to-optics step response 
is governed by the rolloff filter that is 20 dB down at 1 MHz and the 
ideal step function is a good model for the analog system. Since there 
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Fig. 8—SIB vs T'z for different sampling rates. 
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is a filter in the front end of the codec, the step-function model is also 
appropriate for the digital case as long as the frequency characteristic 
of the prefilter dominates that of the camera system. All prefiltering 
that will be seriously studied will be from 6 to 40 dB down at 1 MHz 
and will indeed dominate the camera system. In some instances, re- 
sults will be given for less prefiltering (as in the case of points 0 and 0.1 
in Fig. 7) and then it should be realized that these are somewhat 
idealized cases used only for illustrative purposes. 


3.3 Quantizer 


The addition of the quantizer-decoder combination completes the 
digital system given in Fig. 3. The quantizer is required in order to 
transmit over a finite digital channel. This necessary addition further 
confuses the edge indeterminacy—rise time picture given by the pre- 
filtering and sampling circuits. The quantizer both produces its own 
edge indeterminance and affects average rise times. To see how a quan- 
tizer can produce edge indeterminance, refer to Fig. 9a. Here the step 
response of the prefilter is given by a smooth straight line. The quan- 
tizer used in this example has the standard Phase 0, companded DPCM 
3-bit characteristic described in Ref. 3. With the sampling instants 
given at the bottom of Fig. 9a, the prefilter output could be coded as 
either decoder output No. 1 or No. 2, depending on slight noise varia- 
tions at the beginning of the step. After passing through the receiver 
rolloff filter, the response to decoder outputs No. 1 and No. 2 is given 
in Fig. 9b by rolloff filter outputs No. 1 and No. 2 respectively. The 
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Fig. 9—Quantizer-induced busyness. (a) Prefilter and decoder outputs. (b) Re- 
sulting rolloff filter outputs. 
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ALL STEP RESPONSES FROM: 
1.856-MHz SAMPLER 
PHASE 1—4-BIT QUANTIZER 






0.1 dB 6 dB 12 dB DOWN AT 1 MHz 


Fig. 10—Step response envelopes. 


resulting edge indeterminance (again measured at the 50-percent am- 
plitude point) is due to the quantizer and will henceforth be called 
quantizer-induced busyness (QIB). This quantizer-induced busyness 
is any resulting edge indeterminance over and above that which would 
be produced by just the prefiltering and sampling circuits. The total 
edge indeterminance resulting from a combination of both sampling- 
and quantizer-induced busyness will be referred to as edge busyness 
(EB). In all cases the edge busyness effect is calculated by varying the 
phase of the sampler. This phase difference can result in practice in 
many ways. It occurs spatially for a slanted line. It can result from 
a moving scene, from an unsynchronized codec, or from the residual 
phase jitter in a synchronized codec. It can also be shown that to some 
extent random noise on the video signal can be mapped into an equiva- 
lent sampling phase jitter. 

The envelopes of possible step responses at the output of the receiver 
rolloff filter for the Phase 1 system are given in Fig. 10 for a Gaussian 
prefilter that is 0.1, 6, and 12 dB down at 1 MHz (see the Appendix for 
calculations). This Phase 1 system employs a 1.856-MHz sampler and 
a companded DPCM 4-bit quantizer which is described in detail in 
Table I. Again the prefilter yields a tradeoff between edge busyness 
and rise time. For any combination of sampling rate and quantizer, 
the prefilter yields this same type of tradeoff. To assign numbers to this 
characteristic, Fig. 11 plots edge busyness versus rise time for the 
Phase 1 system with varying amounts of Gaussian prefiltering. Also 
plotted for comparison purposes is the sampling-induced busyness— 
rise time tradeoff already given in Fig. 7. 
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TaBLE [—Puase 1, 4-BiT QUANTIZER CHARACTERISTIC 


Output Level Number Quantized Signal Output* oe) 
1, 2 rae | 
3, 4 + 3 
5, 6 + 7 
7, 8 +15 
9, 10 +23 
11, 12 +31 
13, 14 +39 
15, 16 +47 


* Weights are relative to the smallest quantum step of a 7-bit PCM system. See 
Ref. 4 for additional information. 


The addition of the quantizer complicates the problem even more 
than already discussed. Up to this section the results have been inde- 
pendent of the magnitude of the black-white transition. Fig. 11 results 
from using a 40-percent black-white transition with the Phase 1 4-bit 
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Fig. 11—Edge busyness—tise time tradeoffs. 
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Fig. 12—Edge busyness and rise time vs step size. 


quantizer. Both edge busyness and rise time are plotted versus step 
size in Fig. 12 for the Phase 1 coder using its 6-dB Gaussian prefilter. 
The shape of both curves can be explained using Fig. 18, where both 
the quantizing noise and the step-size-to-quantizing-noise ratio are 
plotted versus step size. Both edge busyness and rise time are inversely 
proportional to the signal-to-noise ratio and are therefore minimized at 
a step size of 30 to 40 percent where the signal-to-noise ratio is maxi- 
mized. For the remainder of this study only three measures of system 
quality will be considered. These three measures are: 


(1) the edge busyness and rise time of the 40-percent step size 
(M4), 
(17) the average edge busyness and average rise time for the 10, 20, 
30, and 40 percent step sizes (M1-4), 
(i201) the average edge busyness and average rise time for the 10 step 
sizes from 10 to 100 percent (M1-10). 
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Fig. 13—Quantizer noise and SNR vs step size. 


Measure M4 has the advantage of simplicity. Edge busyness can 
easily be seen on a transition of this size. This is also the breakpoint 
shown in Fig. 13 before the quantizer’s slope overload characteristic 
takes effect. Measure M1-4 averages over the lower four step sizes. 
Typically 95 to 99 percent of all transitions are below a 40-percent step 
size for a 2-MHz sampler. Measure M1-10 uses the most information 
and disregards all assumed knowledge of observer subjective effects. 
All three measures were used in designing the coders that will be dis- 
cussed in the following sections. Ultimate justification for these mea- 
sures rests with the subjective evaluation of coders with quantizers as 
will be discussed in Section VI. 
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IV. SUBJECTIVE EVALUATION OF THE EDGE BUSYNESS—RISE TIME TRADE- 
OFF 


In the last section it was shown how the variation of the prefilter 
yielded a tradeoff between edge busyness and average rise time for any 
combination of sampling rate and quantizer structure. In this section 
a series of subjective tests will be described that were run to find the 
subjectively optimum edge busyness—rise time tradeoff. For these tests 
only sample-and-hold circuits were used in order to avoid the question 
of which measure to use when quantizers are involved. In Section 
VI a series of subjective tests will be described using coders with 
quantizers. 

The test consisted of asking each of 30 observers to rate various sys- 
tems on the basis of three different source pictures. The first camera 
source was & Mannequin, set up as shown in Fig. 14. This mannequin, 
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Fig. 14—Microhenry. 
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“Microhenry” or ‘‘Mike” for short, is mounted on a platform which 
tilts him in a symmetric, sinusoidal fashion from side to side. He is 
pivoted approximately 14 inches below his chin and was moved a 
maximum of +13 degrees from the vertical at a rate of one cycle every 
4.3 seconds. The second and third camera sources were slides of Karen 
and a graph, as shown in Figs. 15 and 16, respectively. Both of these 
sources were stationary. For all three sources, each of the 30 observers 
was given an AB type test wherein the observer was presented with 36 
different pairs of systems and asked to pick the most pleasing or the 
most readable of each pair (most pleasing for each test given on Mike 
and Karen and most readable for those tests given on the graph). These 
36 pairs include all meaningful comparisons between nine different sys- 
tems. If the pair (A, B) = (system 6, system 9) was given, then pair 
(A, B) = (system 9, system 6) would not be included. Similarly there 
was no AB of any system with itself. Each of these nine systems con- 
sisted only of a prefilter followed by a sample-and-hold circuit. All pre- 
filters had a crispened Gaussian amplitude characteristic which, along 





Fig. 15—Karen. 
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with near-linear phase, resulted in a symmetric step response with 
single 4-percent undershoot and overshoot. Five different prefilters 
were used, being designed to be 3, 6, 10, 20, and 40 dB down at 1 MHz. 
After tuning, these filters were measured to be 3.0, 5.4, 9.6, 18.4, and 
38.1 dB down at 1 MHz with approximately linear phase. The sample- 
and-hold circuit was driven by one of three sampling clocks, each of 
which was synchronized to the camera video. The three sampling 
clocks had synchronized frequencies of 2.015, 1.512, and 1.200 MHz. 
The particular numbering sequence with system characteristics is 
listed in Table II. Both the predicted and the measured performance 
characteristics of these nine systems, as specified on an edge busyness— 
rise time chart, are given in Fig. 17. Each test was given using an 
unmodified Mod 2C Picturephone station set as the receiver, with 
the observer given complete freedom of its brightness control. Process- 
ing the resulting data for each camera source using standard-pair com- 
parison techniques? results in the scale values plotted in Fig. 18. These 
scale values are a psychological measure of relative coder preference; 
the higher the scale value, the more the coder is preferred. The least 
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TasuE IJ]—System CHARACTERISTICS 


System Number Sampling Rate Prefilter* 


1 2.015 MHz 6 
2 1.512 MHz 3 
3 1.512 MHz 10 
4 1.512 MHz 20 
5 1.512 MHz 40 
6 1.200 MHz 3 
7 1.200 MHz 10 
8 1.200 MHz 20 
9 1.200 MHz 40 


* Specified by the number of dB down at 1 MHz. 
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Fig. 17—Edge busyness—rise time tradeoffs for subjective tests. 
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Fig. 18—Psychological scale values for Mike, Karen, and graph. 


preferred coder is arbitrarily assigned the value zero with all other 
coders scaled relative to the worst. The following properties should be 


noted: 


(7) 


(22) 


No intersource numerical comparisons should be made. A scale 
value of 4.13 for the best graphics coder and 3.13 for the best 
coder on Mike only means that the observers were more dis- 
criminating on the graphics source. It does not mean that the 
best graphics coding is better than the best coding of Mike. 
Mike is somewhat of a worst case for aliasing. His shirt was 
chosen for its narrow stripes because undersampling is most 
harmful when appreciable high-frequency content is present. 
This undersampling causes the straight edges to appear crooked 
and the moving subject imparts a distinctive motion to them. 
On Mike most observers found aliasing to be the only degrada- 
tion caused by the different systems and hence the narrower 
prefilters were preferred because they yielded less aliasing at 
the expense of lower picture resolution. A prefilter that is 20 
dB down at 1 MHz was preferred for the 1.5-MHz sampling 
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rate and a prefilter 40 dB down was preferred for the 1.2-MHz 
system. 

(iit) Karen has been a much-used source where both edge busyness 
and resolution are important. The aliasing still shows up as 
crooked stripes in her blouse, but they are motionless for our 
synchronized clocks and therefore not quite as annoying. For 
this source a 20-dB prefilter was preferred for the 1.5-MHz sys- 
tem and a 3-dB prefilter was preferred for the 1.2-MHz system. 

(iv) The graph’s sole property is readability. Here the 10-dB pre- 
filter was preferred for the 1.5-MHz system and the 20-dB 
prefilter was preferred for the 1.2-MHz system. 


Since Mike is a worst case for aliasing, Karen a difficult source for 
both aliasing and resolution, and the graph a test of readability, the 
three sources are considered to be of equal importance, and all 90 re- 
sponses will be lumped together for the final preference scale. Analyz- 
ing these 90 responses by the pair comparison analysis‘ results in the 
scale values shown in Fig. 19. The 2-MHz system is preferred over the 
best 1.5-MHz system which is in turn preferred over the best 1.2-MHz 
system. There is a fairly broad range of optimum prefiltering for the 
1.5-MHz system with the best prefilter 20 dB down at 1 MHz. For 
the 1.2-MHz system, the narrower the prefilter the better, with the 
best choice of those studied being 40 dB down at 1 MHz. 

As shown in Fig. 17, these results indicate that the best subjective 
tradeoff between edge busyness and rise time occurs just below the knee 
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Fig. 19—Psychological scale values resulting from all 90 responses. 
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Fig. 20—Edge busyness—trise time tradeoffs for the Phase 1 quantizer. 


of each tradeoff curve where the incremental reduction in edge busyness 
becomes small compared to the incremental increase in average rise 
time. 


V. OPTIMIZING THE SAMPLING RATE FOR A FIXED TRANSMISSION RATE 


For the Phase 1 4-bit companded DPCM quantizer given in Table 
I, the variation of both prefilter and sampling rate leads to the trade- 
offs between edge busyness and average rise time given in Fig. 20. These 
curves are not unexpected; for any combination of sampling rate and 
quantizer, the prefilter leads to a tradeoff between edge busyness and 
average rise time. The higher the sampling rate (and therefore the 
higher the transmission rate in this case), the closer the knee of 
the tradeoff curve approaches the analog case and therefore the 
better the picture quality. Figure 20 results from using a crispened 
Gaussian prefilter. Many other filter shapes were studied in an attempt 
to find the best edge busyness—rise time tradeoff curves. At the higher 
sampling rates (namely 1.856 MHz and above) the prefilter shape be- 
comes unimportant; here, less prefiltering is required and the receiver 
rolloff filter becomes increasingly dominant. As the sampling rate is 
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reduced, two things are desired in filter shaping: rapid amplitude fall- 
off versus frequency and controlled overshoot in the time domain. The 
rapid falloff is required to control aliasing without removing significant 
video information and controlled overshoot is needed to present a pleas- 
ing picture, one without excessive ringing. The crispened Gaussian 
filter shape with 4-percent overshoots, as specified by eq. (1) with 
K = 0.5292 but with T a variable, yields a good compromise between 
these two factors. 

The results given in Fig. 20 show the effects of different sampling 
rates with a single quantizer. The transmission rate is proportional to 
the sampling rate and the number of bits used to code each sample. 
Therefore, in the extreme case, we are comparing the performance of 
an 8-Mb/s transmission system (2-MHz sampler—4-bit quantizer) 
with that of a 4-Mb/s transmission system (1-MHz sampler—4-bit 
quantizer). The more interesting tradeoff results from considering a 
specific transmission rate. For any given transmission rate, what is the 
best combination of prefilter, sampling rate, and quantizer which 
minimizes both edge busyness and average rise time? A higher sampling 
rate implies fewer quantizing levels which means less SIB but more 
QIB. A lower sampling rate implies more quantizing levels which 
means more SIB but less QIB. The prefilter is used to trade off edge 
busyness against picture resolution to yield the best subjective picture 
quality for each sampling rate—-quantizer combination. 

To be specific, given the following three transmission rates* 


(zt) 6. Mb/s 
(it) 7.5 Mb/s 
(772) 8. — 9. Mb/s 


what is the optimum combination of sampling rate and number of 
bits/sample? 

For the 6.-Mb/s rate, the first three coders listed in Table III are 
considered. In order to simplify the calculations, only the first per- 
formance measure, M4, slightly modified, is calculated. A 40-percent 
video transition is used with the 2.0-MHz system, a 50-percent transi- 
tion with the 1.5-MHz system, and a 60-percent transition with the 


* The 6.-Mb/s rate results from the simplest use of a T2 line. The effective bit rate 
of 7.5 Mb/s results from using 300 cells of storage to expand video information into 
the horizontal sync pulse interval with subsequent transmission over a T2 line. An 
effective bit rate of 8 to 9 Mb/s can be achieved over a T2 line by using about 20,000 
cells of storage (see Ref. 6) to use both the horizontal sync pulse interval and variable 
length coding. This variable length coding takes advantage of the statistical properties 
of the video signal. 
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TABLE IJJ—CopEr CHARACTERISTICS 




















Transmission Rate Sampling Rate Quantizer 
6.0 Mb/s 2.0 MHz 3-bit DPCM 
6.0 Mb/s 1.5 MHz 4-bit DPCM 
6.0 Mb/s 1.2 MHz 5-bit DPCM 
7.5 Mb/s 1.9 MHz 4-bit DPCM 
7.5 Mb/s 1.5 MHz 5-bit DPCM 
7.5 Mb/s 1.2 MHz 6-bit DPCM 
8.3 Mb/s 1.9 MHz 5-bit DPCM 
9.0 Mb/s 1.5 MHz 6-bit DPCM 
8.4 Mb/s 1.2 MHz 7-bit DPCM 











1.2-MHz system. This modification accounts for the longer sampling 
intervals encountered with the lower sampling rates. For each coder 
a wide range of prefilters was evaluated and the optimum prefilter 
chosen using the subjective results given in Section IV. The prefilter 
which results in an edge busyness—rise time tradeoff just below the 
knee of the curve was chosen. The three resulting systems are given by 
the upper curve of Fig. 21. For the effective bit rate of 7.5 Mb/s, the 
middle three coders listed in Table III are considered. Using similar 
optimization procedures, these three coders result in the three edge 
busyness—rise time points given by the middle curve of Fig. 21. For 
the 8. to 9.-Mb/s transmission rate, it is assumed that variable length 
coding will result in an extra bit per sample for each of the sampling 
rates used. Here the three coders listed in the bottom of Table III are 
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Fig. 21—Edge busyness—tise time tradeoffs for fixed transmission rates. 
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considered. These coders result in the three edge busyness—rise time 
points given by the lower curve of Fig. 21. 
Examination of Fig. 21 leads to the following conclusions: 


(t) For intraframe coders operating at an effective bit rate of 6 
Mb/s (points 3b-2., 4b-1.5, and 5b-1.2 on Fig. 21), the edge 
busyness—rise time tradeoff is optimized by using a 1.5-MHz 
sampling rate and a 4-bit quantizer coupled with a Picture- 
phone-type prefilter 20 dB down at 1 MHz. 

(iz) For intraframe coders operating at an effective bit rate of 7.5 
Mb/s (points 4b-1.9, 5b-1.5, and 6b-1.2), the edge busyness— 
rise time tradeoff is optimized by using a sampling rate of 
either 1.5 or 1.9 MHz, with no significant advantage associated 
with either of the two sampling rates. 

(iit) Only after an intraframe coder is operating at a bit rate higher 
than 7.5 Mb/s (points 5b-1.9, 6b-1.5, and 7b-1.2) does it become 
beneficial to sample at 2 MHz. 


VI. SUBJECTIVE EVALUATION OF PICTUREPHONE CODERS 


The main point indicated in Section V was that the optimum sam- 
pling rate is some function of the transmission rate. This point is most 
dramatically illustrated for a transmission rate of 6. Mb/s. Here it is 
predicted that a 4-bit-1.5-MHz coder will outperform the Phase 0 
coder (8-bit-2.0-MHz). In order to check this point subjectively, the 
same 30 observers used in Section IV were given an AB test between 
the Phase 0 coder and a coder consisting of a crispened Gaussian pre- 
filter 20 dB down at 1 MHz, a synchronized 1.512-MHz sampling clock, 
and the Phase 1 4-bit quantizer. Each observer was asked to pick the 
most pleasing coder for Mike and Karen and the most readable coder 
for the graph. A summary of the results for the Phase 0—4-bit—1.5-MHz 
AB comparison is given below: 


Prefer Phase 0 Prefer 4-bit—1.5-MHz 
Mike: D2 8 
Karen: 3 27 
Graph: 7 23 
All: 32 58 


Overall, the 4-bit-1.5-MHz coder was preferred over Phase 0 by a 
margin of nearly 2:1, giving some justification for the two-factor 
tradeoff analysis. 
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In order to find out where the edge busyness—rise time impairments 
lie on an absolute quality scale, a series of unimpaired-impaired type 
comparisons were used to rate six different coders on a five-comment 
impairment scale. For each of the three camera sources, each observer 
made a comparison of all six coders to the unimpaired analog video 
signal and was asked if the impairment added by the coder was: 


(z) not noticeable (weight = 1) 

(iz) just noticeable (weight = 2) 
(77) noticeable but not objectionable (weight = 3) 
(tv) objectionable (weight = 4) 

(v) extremely objectionable (weight = 5). 


Six different coders were evaluated; these were: 


(t) Phase 1 (1.856-MHz, 4-bit DPCM, 6-dB Gaussian pre- 
filter) 
(it) Phase 0 (2.0-MHz, 3-bit DPCM, 6-dB Gaussian prefilter) 
(iit) the 4-bit-1.5-MHz coder used in the previous AB comparison 
with Phase 0 
(iv) a 2.0-MHz sample-and-hold system using a 6-dB crispened 
Gaussian prefilter (System 1 of Section IV) 
(v) a 1.5-MHz sample-and-hold system using a 20-dB crispened 
Gaussian prefilter (System 4 of Section IV) 
(vt) a 1.2-MHz sample-and-hold system using a 40-dB crispened 
Gaussian prefilter (System 9 of Section IV). 


The detailed results of the comment scale ratings in terms of means 
and standard deviations for Mike, Karen, the graph, and all of the re- 
sponses are listed in Table IV. When all 90 responses are lumped to- 
gether, the six coder means are given below (see the next to the last 
row of Table IV). 


Mean Comment Scale Rating 


2.0-MHz sample-and-hold : 1.68 
1.5-MHz sample-and-hold: 2.12 
1.2-MHz sample-and-hold: 3.46 
Phase 1: 2.59 
Phase 0: 3.31 
4-bit-1.5-MHz: 3.30 


Here a lower mean comment scale rating is preferable, indicating a 
lower amount of impairment added by the coder. Note that although 
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TABLE IV—ComMENT SCALE RESULTS 


Coder: 4b-1.5 | Phase 0 | Phase 1 | 2.0 S&H | 1.5 S&H 1.2 S&H 
Prefilter: 20 6 6 6 20 40 dB 
Samp. Rate: 15 2.0 1.9 2.0 1.5 1.2 MHz 
Quantizer: 4 3 4 — — — Bits 

Mike 
Mean: 3.27 2.85 2.50 1.78 2.10 3.05 
Standard 
Deviation: 0.96 0.78 0.86 0.79 0.61 0.78 
Karen 
Mean: 3.23 3.38 2.48 1.28 1.60 3.32 
Standard 
Deviation: 0.63 0.54 0.72 0.50 0.65 0.89 
Graph 
Mean: 3.40 3.68 2.80 1.82 2.67 4.00 
Standard 
Deviation: 0.85 0.66 0.88 0.76 1.06 0.96 
All 90 Responses 
Mean: 3.30 3.31 2.59 1.63 2.12 3.46 
Standard 
Deviation: 0.83 0.75 0.84 0.74 0.91 0.97 


on an AB test the 4-bit-1.5-MHz coder was chosen over the Phase 0 
coder by a margin of almost 2:1, when placed on a quality scale and 
separated in time, they are rated as being equal. In using these ratings 
it should be remembered that the three camera sources were picked to 
give all intraframe coders trouble. If Mike did not have a striped shirt 
or Karen a striped blouse, these tests would have shown very little 
difference between the coders and all would have been rated much more 
leniently. 

In Fig. 22, the six coders used in the above subjective tests are 
plotted on an edge busyness—rise time chart. A 40-percent step size 
is used for both the Phase 0 and Phase 1 coders while a 50-percent 
step size is used for the 4-bit-1.5-MHz coder. All sample-and-hold 
coder coordinates are independent of the step size due to lack of a 
quantizer. With the mean scale values indicated in parentheses, it be- 
comes apparent that the general edge busyness—tise time analysis is 
further substantiated. The closer the coordinates to the analog case 
the better, with almost equal weighting between the edge busyness 
and rise time axes, edge busyness being slightly more undesired than 
reduced rise time. (See Section IV where points just under the knee of 
each curve were chosen.) The two dashed circles with centers on the 
analog coordinates approximate equal-impairment contours with the 
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Fig. 22—Comment scale values as a function of edge busyness and rise time. 


dashed 45-degree line giving equal weighting to both edge busyness 
and increased average rise time. 

Finally it is noted that only station sets using crispened Gaussian 
postfilters which are 20 dB down at 1 MHz have been assumed. Other 
postfilters were tried, including some in the Butterworth family, and 
similar results were obtained. Any postfilter which affects the analog 
system to only a minor degree is not expected to significantly change 
the results presented in this paper. 


VII. SUMMARY 


A description of the two main degradations inherent in present intra- 
frame video coders has been presented. Methods allowing a tradeoff 
between these two degradations, edge busyness and reduced rise time, 
have been described; and optimum coder composition, including pre- 
filtering, sampling rate, and quantizer structure, has been presented as 
a function of transmission rate. A subjective testing program has been 
carried out which substantiates the edge busyness—rise time evaluation 
of video coders and rates a number of intraframe coders on a five- 
comment impairment scale. 
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APPENDIX 
Edge Busyness—Rise Time Program 


Program Inputs: (i) Step Size [A ] 
(ii) Prefilter Step Response [b(é) ] 
(zit) Sampling Rate [f, = 1/T ] 
(iv) Quantizer Structure. 


Program Computations: 
A.1 For each of 20 uniform phasings: 
(¢:,7 = 1, 20) = @ — 1)T/20. 


(3) 


A.1.1 Determine a series of sample-and-hold values for the system 


as given in Fig. 3. 


(Xij, j = 1, 10) = (Ad[¢: + T(j = 1) }). 


(4) 


A.1.2 Transform sample-and-hold series from A.1.1 into a series 
of values that would be transmitted over the digital channel. 


Ya = Q(Xi1) 


(Yi; j = 2; 10) = (2 Ez _ > ¥u|) 


k=1 


(5) 


(6) 


where Q(a) refers to the quantized value of a that can be 
calculated using the quantizer structure. All thresholds are 


located midway between the output levels. 


A.1.3 Digitally add the quantized series from A.1.2 as would be 
done by the decoder to result in the quantized version of the 


sample-and-hold values of A.1.1. 


(Z3, j = 1, 10) = (x Ya) 


k=1 


(7) 


A.1.4 Since the series produced in A.1.3 represents the amplitude 
of a series of rectangular pulses approximating the original 
input video signal, summing a series of appropriately ampli- 
fied and delayed Picturephone receiver step responses as 
given by eq. (2) of the text will result in the video output of 


the Picturephone station set. 


Z(t) = r Yuslt — o: — (k — 1)T]. 


kel 


(8) 
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A.1.5 Use a search procedure on this video output to find the in- 
stants of time when the waveform is equal to 10, 50, and 90 
percent of the input step amplitude. 

A.2.1 Calculate the output 


1 20 
Tr= an a (T iso — Tito) (9) 


as the average 10- to 90-percent rise time. 
A.2.2 Calculate and output 


EB = max (Tio) — min (Tso) (10) 


as the edge busyness measure. 


The above procedure can be used to calculate the combined effects 


of SIB and QIB or the effect of SIB only, since replacing the quantizer 
with an identity equation leaves only the effects of SIB to be calcu- 
lated. It is not possible to calculate the effects of QIB only, since 
presampling is required for the quantizer. 


For all plots given in this paper, eq. (8) can be placed inside a DO 


loop to calculate the output waveform for a wide range of times. 
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Synthesis of Multiple-Feedback 
Active Filters 


By G. SZENTIRMAI 
(Manuscript received October 3, 1972) 


A synthesis technique is developed for the active RC realization of 
transfer functions that have all their transmission zeros on the tmaginary 
(jw) axis. The method leads to a realization using biquadratic blocks in a 
multtple-feedback arrangement that 1s the generalization of the structure 
obtained from the passive, double-terminated reactive equivalents. The 
realization combines the easy independent tuning properties of the cas- 
cade, with the low-sensitivity characteristics of passive ladders. 


I. INTRODUCTION 


In the last two decades, the question of how to realize a prescribed 
rational transfer function 
N(s) = Vout 
T(s) = = — (1) 
D(s) Vin 
by RC active structures was the subject of more than a thousand 
learned papers. The consensus at the present time seems to be as 
follows: 





(t) The active element to be used is an operational amplifier. 
(zz) Subnetworks (building blocks) realizing biquadratic transfer 
functions 


no + nis + njs? —N;(s) 
di + dis + dis? _D,(s) 








Ti(s) = (2) 


are constructed as intermediate steps. 
(i7t) Finally, the overall transfer function is realized as a cascade 
connection of these leading to: 


T(s) = Tl T',(s). 


527 
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Our objective in this paper is to challenge step (727) above and in- 
troduce an alternative synthesis method valid for a large class of 
transfer functions. This is done by adding one or more feedback loops 
to the cascade configuration. The effect of this will be that, while the 
transmission zeros of the overall system will remain the concatena- 
tion of those of the individual biquadratic blocks, this will not be 
true for the poles any more. 

We will assume (7) and (zz) to be correct without, however, worrying 
about the details of the specific configuration to be used in step (72) 
above; that is to say, our most elementary building blocks will be 
“black boxes” realizing transfer functions of the form given by (2). 
Next we select a particular structure for our investigation and motivate 
this selection on the basis of prior work. This is followed by the de- 
velopment of a synthesis method for the selected structure that 
simultaneously proves the generality of it. 

Finally we will illustrate the method and demonstrate its advantages 
by an example. 

To circumscribe the class of problems that we will consider, note 
that T(s) in eq. (1) is a real, rational function of the complex fre- 
quency variable s, and therefore both D(s) and N(s) are real poly- 
nomials with the further restrictions: 


(i) D(s) is a strict Hurwitz polynomial, i.e., its zeros are in the 
open left half of the s plane, 
(17) the degree of N(s) is not greater than that of D(s), in notation: 


6N S 6D. 
We will need the additional restriction: 


(iz) all the zeros of N(s) are on the imaginary (jw) axis of the 
s-plane. As a consequence, N (s) is either pure even or pure odd. 


This restriction, while quite serious, still leaves a very large group 
of useful functions to be considered, especially if the reader notes that 
the same restriction applies to functions realized by passive lossless 
ladders devoid of (magnetic) coupling. 

The reason we are searching for new configurations can be traced 
back to a paper by Orchard,! where he has shown by a very simple 
physical argument that double-terminated lossless passive structures 
have very low sensitivity to changes in component values inside the 
passband where the loss is near zero. Furthermore, it is also known 
that ladders are good for maintaining high out-of-band suppression. 
While Orchard’s argument is not completely valid for voltage transfer 
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Fig. 1—Passive ladder. 


functions, and therefore cannot be transferred to active realizations, 
improved performance has been obtained by simulating passive 
ladders by active structures. 

The first results of this kind were those of Girling and Good? and 
Adams.’ Girling and Good recognized that, say, a ladder low-pass filter 
of the form shown in Fig. 1 is described by the equations (also see 
Ref. 4): 


1 
Ty = —— (Vin — V2), 
: Bat Le ) 
1 
Ve = — (1 — Is), (3) 
Cos 


1 
I; = Les (V2 — Vout), 


and 
Vout = eee I 
out Cus ae Gy 3) 

and, as such, are also the describing equations for the active RC-struc- 
ture of Fig. 2. When it came to more complex filters, the structures 
advocated by Girling and Good become considerably more complex 
with nonconstant feedback and other undesirable features. 

Adams,’ working independently, went one step further and realized 
that the standard low-pass-to-bandpass transformation applied to 
both Figs. 1 and 2 would result in the very convenient realization 





Fig. 2—Simulated passive ladder. 
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Fig. 3—Multiple feedback bandpass. 


of Fig. 3 for bandpass filters. Note that all blocks have biquadratic 
transfer functions 7';, and further note that 7, and 7, (the first and 
last ones) have poles of finite Q, but all intermediate blocks have poles 
on the imaginary axis of the s-plane. These blocks are therefore, 
strictly speaking, unstable standing alone, but the overall structure 
is still stable. This becomes less of a surprise if we note that ideal 
reactive elements (L’s and C’s) used in passive synthesis are also, 
strictly speaking, unstable. The only difference is that natural dis- 
sipation will always push passive elements toward the stable side of 
the s-plane while marginally stable active building blocks may hover 
over either side of the jw axis. 

Adams’ numerical results were sufficiently encouraging, showing 
substantial reduction in sensitivities to start him and others on the 
road searching for similar realizations for more general filter functions. 

There are many ways to handle this problem, brute-force numerical 
matching, flowgraph manipulation, and matrix operations on the 
state-variable equations being some of them.® All of these were used 
by researchers at one time or another, with mixed results. They were 
unable to explain the origin and uses of multiple solutions; or if they 
were able to, they led to nonminimal realizations or realizations with 
much more complex feedback and feedforward paths.® 





Fig. 4—Passive bandpass filter. 
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This author also joined in the search and, to his delight, found a 
group of bandpass functions where the approach of Figs. 1 and 2 was 
still working. Consider the passive bandpass filter shown in Jig. 4, 
which is a so-called minimal-inductance realization of an even degree, 
antimetrical filter. This network is described by the equations: 


fy-= YiWVia— Vs) 
Vo= Z2(1i ie Is) 


(4) 
I; = Y3(V2 ane Vent) 
Vout = Zals 
or, with a slight modification: 
I Y 
= = = (Vin — V2) 
s 8 
I I 
Ve a sZ2 (= = =) 
8 s (4a) 
Iz; Ys 
oS ae (V2 — Vout) 
s S 
Is 
Vout = sZ4— 
S 


and this is identical to the equations describing the structure of Fig. 3, 
if we identify: 
Y, Y3 
Ty =— } T3 cy 
s s (5) 
T's = sZa; T's = sZ4. 


Note that all 7; are again biquadratic, 7, and 7 have poles inside 
the left-half s-plane, while the poles of T, and 7’; are on the imaginary 
axis. 

This led us to expect the general structure of the form of Fig. 3 
to be in some ways canonical. Note that the restriction of all feedback 
coefficients to —1 does not restrict generality. For instance, the center 
loop feedback coefficient may be changed to —k if we simultaneously 
multiply 7; by k and divide T, by k, with no change in the overall 
transfer function. This procedure, in fact, can be used to adjust 
voltage levels inside the filter to optimize dynamic range at a later 
stage in the design process. 
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II. SYNTHESIS PROCEDURE 


Consider now the general structure shown in Fig. 5 with the following 
assumptions: 


(z) All 7; transfer functions are real, biquadratic functions of s. 
If the overall degree is odd, 7; or T', will become bilinear. 
(iz) All internal block transfer functions 


T2, Ts, aes T m1 


are even functions of s, i.e., all their singularities are located 
symmetrically either on the imaginary or the real axes. 

(zit) The poles of 7; and T,, are inside the left-half s-plane, and 
their numerators are either even or odd. 


The synthesis problem can now be formulated as follows. Given a 
suitably restricted overall rational transfer function of the form of 
eq. (1), find the biquadratic transfer functions 


= N.(s) 
7 Dj(s) 


satisfying (2) through (777) above, such that these in the structure of 
Fig. 5 realize T(s). 

The outline of the proposed solution of this problem is as follows. 
First we recognize that our structure is equivalent to a (reciprocal) 
ladder network, if we restrict all the feedback coefficients to be —1. 
This restriction will be removed later in Appendix B. The synthesis 
problem can therefore be solved if we somehow derive a complete 
(impedance, admittance, or any other) set of parameters for this 
ladder. Once a set of parameters is obtained, the actual synthesis can 
follow along lines very similar to the well-known passive filter synthe- 
sis method. 

Let us use our ladder — multiple-feedback analogy backwards now 
and construct the ladder network of Fig. 6a (assuming that m is even), 





T :(s) 7=1,2,---,m (6) 





Fig. 5—General multiple-feedback structure. 
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A A A 
Y,=Ty Y3=T3 Yona = Tm 





Fig. 6—(a) Equivalent pseudo-ladder for m = even. (b) Equivalent pseudo-ladder 
for m = odd. 


where the branches need not be realizable passive immittances. How- 
ever, the voltage transfer functions of the two structures are clearly 
identical and must be equal to T = N(s)/D(s). Note first that 
obviously : 


N(s) = tl N(s). (7) 


i=l 


The Z and Y matrices of our pseudo-ladder must be of the forms: 





1 TD N 
Z= (8a) 
NiP LN NiN»R 
1 [NiN,R —N 
ae EE (8b) 
Nr»Ql —N D 


where P, Q and R are three unknown polynomials. These follow from 
Fig. 6a and the assumptions (2) through (77) above. Considering 
degrees and parities of the polynomials in question, we see that, as- 
suming N, and N,, to both be even: 


5D = 2m 6P = 2(m — 1) 
5N = 2m and pure even 6Q = 2(m — 1) (9) 
5R = 2(m — 2) and pure even 
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where 6X, as before, denotes the degree of the polynomial X. That R(s) 
must be even follows from the fact that, with a short-circuited output 
I, = You N N 





qT; Vo=0 Yu NiN»AR 7 R 


and since this parameter is independent of both Y,; and Zm, it must 
be pure even. Here we used the notation: 


m—1 


N= TI Ni. (10) 


t=2 


Since ZY = 1, the 2 X 2 unit matrix, we get the determinantal 
relationship first: ; 
DR — N,N,N? = PQ. (11) 


Note that N; and N,,, must be known at this stage. The significance 
of this is explored further in Appendix B. 

In order to proceed further, we must consider the actual synthesis 
procedure. Consider, for instance, the parameter 


Nik 
Q 
We can get Y, by extracting the partial fraction 


ai + Bis 
N; 





Yu = 


plus a constant from 1/yi: such that the remainder has a factor 
N, in its numerator. More importantly, the remainder must be a pure 
even function of s, since it 1s now independent of both Y; and Zn. 

Separating Q(s) into even and odd parts and performing the opera- 
tions outlined, we see that Q(s) must be of the form: 


ki _- 
Q(s) = A(s) + S sR(s) (12) 


where 6A = 2(m — 1) and pure even, and k, = 26; is a constant. A 
similar argument applied to z22 shows that P(s) must be of the form: 


P(s) = B(s) + = sh) (13) 


where 6B = 2(m — 1) and pure even, and k,, is a constant. Let us 
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now separate D(s) into even H2(s) and odd s£,(s) parts: 
D(s) = sEy(s) + E2(s). (14) 
Substitute (12) through (14) into the determinantal equation (11) and 
separate even and odd parts: 
2h, =k»A +hiB (15) 


‘ k? 
E2R — NiN,N? = AB + Faas (16) 


where we introduced kyk,, = k?. 
Considering eq. (16) a quadratic in the unknown polynomial R(s), 
this will have a polynomial solution if and only if the discriminant: 


E3 — k’s(AB + NiN»N?) © @? (17) 
is a full square for some even polynomial G. In such a case, the solu- 
tion is given by: 


2 
Tr, ees, 


and, in order to get a polynomial, we must select the negative sign, 
since from (17) we see that the constant terms of /, and G will be 
identical and hence cancel from the difference: 


_ 2B — G) 


k2s2 


R (18) 


Let us next express, say, B from eq. (15) and substitute it into (17): 
E3 — k?s?NiN,N? — 2kms*HiA + k3, s°A? = G?. (19) 


This is again a quadratic in A, having a polynomial solution if and 
only if : : 
s?H? — £3 + k’s?NiN,,N? + G? = s*H? (20) 


for some even polynomial H. Note that, by eq. (17), the left side of 
(20) must have a double zero at s = 0; hence the factor s? on the right. 
The solution to eq. (19) is then of the form: 





FitdH 
A=— (21) 
Kin 
and consequently : 
Ei, +H 
i . (22) 





ky 
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Now we have exchanged our three unknown polynomials R, A, and B 
for two new ones, G and H, both even and 

6G = 2m 

6H = 2(m — 1). 
Let us however write eq. (20) in the form: 

(E3 — s*E?) — k’s?N.N,,N? = (G? — s*H?). (23) 
Now introduce the notation: 
F(s) = G(s) + sH(s) 

and use (14) to write eq. (23) in the form: 


D(s)D(—s) — P?8:N.N,,N? = F(s)F(—s) (24) 





which is our design equation. Once we pick N,; and Np, the left side is 
known save for the constant k?, which can be selected arbitrarily as 
long as the left side does not have pure imaginary roots of odd mul- 
tiplicity. This condition is clearly satisfied for k? = 0 and hence, by 
continuity, there must be a finite range 


2 


OSk Skinax 


IA 


for which it is satisfied. 
Our synthesis process is now straightforward and is as follows: 


(z) Select Ni and Np. 
(21) Select k? arbitrarily but such that 


D(s)D(—s) — k?s’N,N,,N? 


has no imaginary roots of odd multiplicity and factor this 
polynomial into 
F(s)F(—s). 


Note that F(s) need not be a Hurwitz polynomial. 
(iz) Generate G and H from: 


F(s) = G(s) + sH(s) 
and EF, and E, from: 
D(s) = E2(s) + sE,(s) 


where £, H,, G and H are all even polynomials. 


(wv) 


(v 


nN 
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Form A, B and R from eqs. (21), (22) and (18) respectively 
and P and Q from eqs. (18) and (12), giving: 


(EH. — G) + s(#i + A) = D(s) — F(+s) 


P= (25) 
kis kis 
on (By — G) + s(B1 H) — D(s) — F(¥s) 
Q= her ee es ee, (26) 
kins kn8 


Select the sequence of the remaining numerators N;(s), 7 = 2, 
3, --:, m — 1 and synthesize the network by a technique very 
similar to the standard “zero shifting’ technique’ used in 
passive reactive ladder network synthesis. The difference is 
merely that here we deal with even rather than odd rational 
fractions. 


Clearly, this process has several problems. 


(a) 


(b) 


(c) 


The above is only one of many cases depending on the parity 
of m, the parity of D, N,, and N,. All these cases have been 
analyzed, and the general design equations are summarized in 
Appendix A. 

No guarantee is available that the resulting 7;(s) functions 
will have all positive coefficients. This is similar to the passive 
ladder form of realization, where all positive element values 
cannot, in general, be guaranteed either. However, positive 
coefficients for 7;(s) are not necessary for realization. See 
Appendix D for further comments on this matter. 

We must make the more-or-less arbitrary choices in steps (2) 
and (iz) above. Our first concern is to minimize overall sen- 
sitivity. We will consider this problem in Appendix B after we 
have discussed the calculation of the sensitivity. 


III. CALCULATION OF 7’; TRANSFER FUNCTIONS 


Step (v) above concerns the final step, that is, the calculation of the 
coefficients of the individual 7; transfer functions. The method is a 
slightly modified version of the reactive ladder synthesis, but for 
completeness we describe it here in some detail. 

First we note that we must select an arbitrary (positive) constant 
k, and then k,, is determined from kik, = k?. This constant is im- 
material, since we may later multiply all odd indexed 7; by a constant 
and divide all even indexed 7; by the same constant. This procedure 
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will leave the overall transfer function unchanged if there are an even 
number of 7’; blocks, or will multiply it by the same constant if there 
are an odd number of blocks in the network. 

For computational simplicity, let us start by using y1:, which in the 
case of the network shown in I’ig. 6a does not depend on Zn, hence 
will only suffice for the calculation of 71, Ts, ---, T’m—1. The last block 
will therefore have to be calculated from another function, say, 229. 
From eqs. (8b), (18) and (26) we get: 


Nik 2Ni(E2—-G@) 
Q — kys[D(s) — F(+s)] 


Inverting this, we immediately recognize that it can be written in 
the form: 

(3) 1 kis[ Ee + shy =. G xt sH | kis ky s*( By + H) 
4 —S=S=S— C—O aC — 
Yu 2Ni(E2 — G) 2Ni 2 Ni(E = G) 


(27) 


Yu = 


(28) 


Both of these terms are known and the second is pure even; hence 
from this point on (that is, nearly from the beginning) we are dealing 
with pure even functions simplifying the work considerably. Denoting 


_ ki s°(Hi & H) 
2 Ni(F. — G) 
we first note that this still has a pole pair at s = + juwi, the zero of 


Ni, because of the factor N; in its denominator. We remove this by 
calculating : 


(it) 2r1 (29) 


ky st(ky is oo H) 


111 = Nz, = 30 
( a, oa 8=jW, 2 Ee —G s?=—w} ( 
and calculate: 
ay 
er2 = 2rl — Ny (31) 


that will not have the factor N; in its denominator any more. Next we 
select the transmission zero + jwe for the next block and calculate: 


(iv) bi = Zeal aces (82) 


which leaves a remainder: 
er3 = &r2 — Bi (33) 


that will have the factor N. in its numerator. At this stage we com- 
pleted a full cycle and calculate the transfer function of the first block 
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as: 


ky 
—_S§ 
(v) eee ae at (34) 
Vv —_—>S  =— + _ = — 
PN Ne 


and have a remainder z,;3 that is two degrees lower. We may invert 
this now and go back to step (77) (with an admittance function this 
time) and repeat this loop as many times as required. Equation (34), 
of course, will contain only the a;/N; and 6; terms giving us the re- 
quired even 7’; transfer functions except when 7 = 1. 

Finally, to obtain the last block as well as to check on the ac- 
curacy of computations, we repeat the same process from the output 
end using 22. The resulting two structures should be identical save 
for a constant multiplier, that is to say, denoting the 7; transfer 
functions obtained from z22 by T;, to distinguish them from those (7) 
obtained from y11, we must have: 


{ee 

T, = CT» 

,-—T 

ae 6. 3 

Te=aHCTs (35) 
T it 

a 6. m—1 

CT. 


T, and 7, have no pairs since y11 and 22 are independent of 7’ and 
T, respectively. The degree to which this set of equations is satisfied 
with a constant C is the numerical accuracy maintained throughout the 
computation. The procedure varies slightly when N(s) is odd and/or 
the degree of D(s) is odd, but only in calculating the first (and last) 
block. These variations are self-evident and hence will not be dealt 
with here. 

Clearly, the internal sequence of zeros (those excluding N; and Nn) 
is arbitrary at this stage. 

Finally, one may note that to ease the numerical accuracy problem, 
the commonly used transform variable®? or the product method” may 
also be used. 
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IV. EXAMPLE 


Many numerical examples have been calculated by the afore- 
mentioned method; here we will illustrate it with one. The example is 
a bandpass filter with passband from w,4 = 0.8 to ws = 1.25 with a 
0.5 dB passband ripple. The filter contains four pairs of transmission 
zeros at: 


Zi,2 = + j0.25 

Z3,4 = + j0.50 

Zs5,6 = + j2.0 

21,3 = + 74.0. 

The corresponding transfer function poles are at: 

Py. = — 0.028950107 + j0.79624226 
P3,4 = — 0.087386703 + j0.90192127 
Ps,6 = — 0.10642659 + 71.0984326 
P73 = — 0.045602221 + 71.2542411. 


Following our procedure, the resulting structure contains four bi- 
quadratic blocks with transfer functions of the form of eq. (2); the 
coefficients are tabulated in Table I. 

For the synthesis, we have selected (see Appendix B): 


N, = (s? + (0.25)?) 
Nm = Na = (s? + (4.)?) 
and 
k? = 3.3 * 10-6 


which is slightly below k2,,. Next, we selected ki = k and kn =k 
and performed the synthesis in both directions. These had a ratio 


C = 8.00488 X 10° 


that was constant to the indicated six decimal digits through all the 
coefficients. Rescaling the coefficients by VC (that is to say, selecting 
k, = k/VC and k, = VCk) resulted in identical values up to six 
decimal digits in the two syntheses, and these are the values tabulated 
in Table I. The computed performance of this structure is shown in 
Fig. 7. All computations were done in double precision arithmetic 
(~ 16 decimal digits), and the indication is that, for higher order cases, 
either the transformed variable method" or, preferably, the product 
method” will have to be used to avoid numerical accuracy problems. 
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TABLE I—FINAL COEFFICIENTS OF EXAMPLE 


t 1 2 3 4 
No 0.0625 4.0 0.25 16.0 
m1 0 0 0 0 
Ne 1.0 1.0 1.0 1.0 

do 7.98786 3.08873 9.31590 11.3657 
dy 2.56983 0 0 2.56983 
de 7.74668 3.02573 9.50989 11.6790 


4.1 Practical Results 


The filter, scaled to 1-kHz center frequency, was constructed in the 
laboratory, using the three-operational amplifier biquad realization” 
for the second-order building blocks. The total structure (see Fig. 8) 
needed two additional phase inverters to get the correct signs for the 
feedback loops. The component values are shown in Table IJ. The 
measured loss (save for a 3-dB flat difference) is also shown on Fig. 7. 
The agreement is the kind or better than the kind one usually expects 
in a passive realization, with no tuning. 


—— COMPUTED 
Oo MEASURED 


DECIBELS 





0 0.5 1.0 1.5 2.0 2.5 
FREQUENCY IN kHz 


Fig. 7—Computed performance of the filter of example. 
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100 KQ. 


R3 


100 KD. 


Fig. 8—(a) Basic second-order block. (b) Complete filter. 


TABLE II—ELEMENT VALUES OF THE REALIZATION OF IIa. 9 


48.1 ko 
15.4 ko 
16.7 ko 
392. ko 
1.97 Ma 
77.7 ka 
9.53 koa 
10. ka 


0.01 uF 
0.01 uF 


15.8 ko 
15.8 ka 


12.1 kQ 


Sections 


6.42 ko 
6.42 kQ 


14.5 kQ 
3.29 ka 
3.16 ka 

169. ko 
2.34 kQ 

17.2 kQ 
1.47 ka 
1.47 ka 


0.05 uF 
0.05 uF 
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V. SENSITIVITY 


There are many ways to analyze the sensitivity of different realiza- 
tions, and the final component sensitivities will clearly depend on the 
particular realization used for the individual biquadratic blocks. In 
order to avoid this last step and provide some meaningful insight into 
the sensitivity properties of the multiple-feedback structure, we will 
compare it with a straightforward cascade realization. 

If several biquadratic sections with transfer functions 7';(s) are 
connected in cascade, the overall function is given by: 


T(s) = IJ Ti(s) (36) 
t=1 
and consequently the sensitivities are simple: 
T T; oT 
Sr; = ——— =1 for all:z. (37) 
T oT; 


Clearly in our multiple-feedback structure, the situation is somewhat 
more complex. However, by the use of continuants," we can derive 
closed-form expressions. 

A continuant K,(X1, +--+, Xn-2, Xn_1, Xn) is defined by the re- 
cursion formula: 
Kn(X1, yey S X pict; Xn) = XnKni(X1, aa Kgs X n-1) 

ae Kn-2 (Xi, eer) X n—2) (38) 

and the two starting values: 


Ki(X1) = Xi. oe 


With these definitions, the overall transfer function of the multiple- 
feedback structure (with all feedback coefficients equal to —1) is 


given by: 
1 1 1 1 
—= Kn, —,—,°:',—}> (40) 


Differentiating and rearranging, using some of the properties of 
continuants, we get the following expression for the sensitivity : 


ToT T Kn i 
aie - (41) 
Tor: T, 1 Lex 

T; 











T 
Sr; —, 
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where 





- é 1 ‘ 1 L) 
Roepe < Pipes pa > Sa 
re eS ae Se ee, AD) 


1 1 1 1 
Ce) See 
T1 Dey T isa Tm 


From this expression, we see that near a transmission zero, where 
T,; = 0, the sensitivity is about unity, that is, it is the same as for the 
cascade case. However, wherever 7; > 1, which will occur somewhere 
in the passband, the corresponding sensitivity becomes very small. 
No other feature is obvious from this equation, hence a specific case 
will be used to illustrate the numbers involved. 

One additional problem arises in the case of the multiple-feedback 
structure, namely, the feedback coefficients can also vary from their 
nominal (—1) value. This can be taken into account by substituting 
a change 





—-1l->- (+8) 
in the kth feedback path by the changes: 


T k4-2441 <= (1 + €) T rpoiyt = 0, 1, coe 
and 
Trpoi— (1 + ©) TT eye: 4=1,2,--:. 


This way we do not need to develop additional formulas for the sen- 
sitivities of the feedback coefficients. 

The resulting sensitivity, neglecting the possible multiplier on the 
overall transfer function, can therefore be derived as follows: A 
change AF;,/F;, in the kth feedback coefficient has the same effect, 
AT, on the overall transfer function as the changes: 




















AT n+: F AF, ; 
——" = (—1)#1 4=1,2,---,m—k. (43) 
Peri Fy, 
In other words: 
AT Ss AT — ( ph la meal (44) 
7 iat OL E45 a i=l OT n 45 e F, 
or 
m—k : 
AT = AM ~ (-) Piss OL 
T Fr in LOT 3 
AF, m—k i 
=—— > (-1)""'Sh, (45) 


F, i=1 
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and finally: 


t+1 T 


m—k 
Sie (CED Shy (46) 
i=1 


Consequently, the sensitivity to a feedback coefficient can be easily 
calculated in terms of the already calculated sensitivities to the second- 
order block transfer functions. 

In order to simplify the problem of comparing the multiple-feedback 
structure with, say, the cascade realization, we have adopted the 
following procedure. 

We calculated the sensitivities of the overall transfer function T 
with respect to all the coefficients c; of the individual biquadratic 
blocks T;. This was done for both the cascade and the leapfrog- 
feedback cases, and in the latter the parameters c; also included the 
feedback coefficients. 

Next we calculated the quantity : 


o = > (Re ST) (47) 


as a function of frequency. This quantity would give a measure of the 
spread of the loss if all coefficients are assumed to be statistically in- 
dependent variables with the same standard deviation. 
Finally, the quantity 
oo 
R= 10 logio a (48) 


MF 


is calculated where the subscript MF indicates multiple-feedback, 
while the subscript C refers to the cascade configuration. R is there- 
fore a measure of the improvement of the sensitivity of the multiple- 
feedback structure over that of the cascade realization. 

This quantity R is plotted as a function of frequency for the ex- 
ample of Section IV above and is shown in Fig. 9. The curve is ex- 
tremely interesting, and it shows a slight deterioration in this measure 
of sensitivity in the stop-band, but a spectacular improvement in the 
passband. The worsening in the stop-bands is minor, but the pass- 
band improvement is substantially more than an order of magnitude. 
This behavior is very desirable since a 1-dB spread in the stop-band 
is almost always immaterial, while the same spread in the passband 
loss can be disastrous. 

In order to provide another comparison, a Monte Carlo yield study 
was made of this filter under the following assumptions. 
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20 


10 LOG, 9 


0 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 
FREQUENCY IN kHz 


Fig. 9—Sensitivity improvement. 


(1) Passband is acceptable if loss varies less than 1 dB from 800 Hz 


to 1250 Hz. 
(iz) Stop-band is acceptable if loss is at least 60 dB from 0 to 500 


Hz and from 1310 Hz up. 
(iit) The operational amplifiers are close to ideal. 
(iv) All passive components have the same tolerance with a flat 


distribution. 
The resulting yields are tabulated in Table III for various tolerances. 


Tor comparison, the same calculations were performed for a cascade 
realization of the same filter, under the same assumptions. The results 


speak for themselves. 


TABLE IJI—ComPpaRISON oF YIELDS OF EXAMPLE 
AND EQUIVALENT CASCADE 


Yield % 
Tolerance MF Filter Cascade 
1.00% 44 ae 
0.50% 78.5 25 


0.25% 100 67.5 
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Clearly, one cannot draw conclusions on the basis of a single ex- 
ample, but a general proof of the low sensitivity of the multi-feedback 
structure is lacking. Note also that no attempt has been made in the 
design of this example to minimize sensitivity apart from the crude 
heuristic arguments indicated in Appendix B. 


VI. CONCLUSIONS 


We have demonstrated that the multiple-feedback structure of 
Tig. 5 is a general one for transfer functions with pure imaginary 
zeros, and we developed a general synthesis procedure for this struc- 
ture. The advantages of this realization are the low-sensitivity prop- 
erties of passive ladders combined with the individually tunable 
transmission zeros of the cascade structure. The resulting configura- 
tion is minimal (no pole-zero cancellation occurs) and hence also 
stable. 

The low-sensitivity properties of the realization have only been 
demonstrated through examples; a general proof of it is still to be 
found. Furthermore, apart from some heuristic arguments, no general 
guidelines are available for selecting the one continuously variable 
free parameter or to help us at the two additional places where dis- 
crete choices are to be made. 

Apart from, conceivably, further optimizing the sensitivity prop- 
erties of the structure, one could also use these choices, as well as 
additional scaling of the individual biquadratic transfer functions and 
the feedback coefficients, to optimize the dynamic range of the filter. 
All of these questions merit further investigations. 
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APPENDIX A 


Summary of General Results 


The results given in Section II of this paper were derived for a 
special case. The general results are as follows: 
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Given 
N(s) 


Tea) 


where the degree of D is n, and denoting: 


n-+ 1 
m= | (number of blocks) 


i.e., the integer part of 1/2(n+ 1) and 


Ms) = > Ni(s) (49) 
N(s) = Ni(s)Nn(s)M(s) (50) 


where all N;(s) are, at most, second order, and N,(s) and N»(s) are 
either even or odd, while the others are strictly even. Now we define: 





Ni,m(s) if Ni,m 1s even 
t= 2 51 
Eiger ie (1 hp ) if Nim is odd. on 
Ww? 
l,m 

The defining equation for 

F(s) = G(s) + sH(s) (52) 
is: 


D(s)D(—s) — (—1)k?s?V1(s) Vn(s)N?2(s) = F(s)F(—s) |. (53) 


The resulting pseudo-ladder is of the form shown in Fig. 6a if ‘‘m’’ is 
even, while if ‘‘m”’ is odd, it is as Fig. 6b, where: 











Nifs) 
T;(s) = 4= 1, 2, ,m 
D;(s) 
The Z and Y matrices of these networks are of the form: 
1 TD N 
“ele tel 
NiP N NiNaAR 2 
1 [NiN,R —N 
Y= | | (55) 
NwQ —N D 
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for ‘“m’”’ even and 


1 NnQ N 
LS ee (56) 
NWrk Ll N NiP 
1 [NiP —N 
Y= =| (57) 
D —N NnQ 


where ‘‘m” is odd. In these equations, the still undefined quantities 
are the polynomial R = R(s) that is always pure even, and the poly- 
nomials P = P(s) and Q(s) that are given by 

a 7 D(s) & F(+s) 


gis)} ~~ kay 2°) x08) 


and 
R(s) = Ev ee + F(s) ] (59) 
k?s2W4(s) W m(s) 
where Ev{ } means the even part of { }, 


Ni,m(s) if Ni,m(s) is even 


— 2 
Wim(s) = (1 + —) if Ni,m(s) is odd om 
(63) 


l,m 


and 
kikm = k, 


Here as well as above in the definition of V1,m(s) the factors: 


s? 
(1 a? 


1,m 





) = Ev D,,,(s). 


Hence, they are unknown initially. Note also that if n is odd, then one 
of the D; and D» is only linear in s, i.e., the corresponding w; or w*, > ©. 

For that case where w; and/or wm are present, an iterative procedure 
must be employed, since these factors must be such that P(s), Q(s) 
and R(s) are all polynomials of the correct degree. 

This condition is satisfied if (1 + s?/w%) and (1 + s?/w?) are factors 
of (#1 + H) and (£, — #) respectively, and thus can be obtained 
iteratively as follows. One picks the factors (1 + s?/w;) and (1 + s?/w”,) 
arbitrarily, but their zeros should be in or near the passband. Then 
one performs the factorization of the polynomial F(s) and calculates 
the zeros of both (#, + H) and (£, — H). At this stage, we can re- 
place our arbitrary factors by the nearest factors of (#, + H) and 
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(E, — H), selecting one factor from each. A repetition of this process, 
or a modification of it, will converge to the desired solution. 

If there is only one unknown factor involved, the iteration is some- 
what simpler and any of the factors of (H; + H) can be used. One 
must further show that, once this iteration converged, one can always 
select sign combinations in the expressions for P(s), Q(s) and R(s) 
such that the resulting expressions reduce to simple polynomials of 
the correct degrees. This can be shown directly in the limiting case 
k? = 0 (when, of course, no iteration is needed), and the existence of a 
solution for a (finite) range of nonzero k? > 0 values can then be 
inferred from continuity again. In order to keep the length of this 
paper within bounds, the details of this step are left for the reader. 
Finally, note that not all sign combinations shown in eqs. (58) and 
(59) are allowed, but at least one will always work. 


APPENDIX B 
Selection of the Parameter k? as Well as the End Factors Ni and Nm 


In Sections II and III we pointed out a number of arbitrary choices 
the designer must make. The following comments may help in some 
of these decisions. 

In order that eq. (53) be factorable in the form of F'(s)F(—s), with 
F(s) a real polynomial, it is necessary and sufficient that the left side 
have no imaginary root of odd multiplicity. Since this is clearly true 
for the case k? = 0 and since the roots will be continuous functions of 
k?, the condition will be satisfied for a range of values 


OS Sx: 


The upper limit k2,, can be calculated as follows: Since D(s) cannot 
have a pure imaginary root, F'(s)F(—s) will have no imaginary root of 
any multiplicity as long as 

F(s)F(—s) je ee Vin(s) N?(s) 
D(s)D(—s) i: D(s)D(—s) 





i” >0 (61) 


s=jw 
or as long as 


, P= Dts?Vi(s) Vin(s) W?(s) 
‘ | D(s)D(—s) io =a 


and therefore this condition will be satisfied as long as 
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wo 


In the limiting case k? = ka there will be at least one s = jw value 
where F'(s)F(—s) will be zero, but all of these will necessarily be of 
even multiplicity. Moreover, since this limiting case in passive filters 
corresponds to the maximum power transfer from source to load, we 
conjecture that it will also be optimum in our active case in the sense 
of providing us with the least sensitive results. 

Going one step further, we note that in the passive case all the pure 
imaginary zeros of F(s) fall inside the passband. In order to achieve 
this in the active case, as well as to make as many zeros of F(s) ap- 
proach the imaginary axis as close as possible, the factors Vi(s) and 
V m(s) should be selected such that 


(—1)™s?Vi(s) Vin(s) 


ss (63) 
Ni(s)Ni(—s8)Nn(s)Nm(—S) | sajo 


is as close to a positive constant in the passband as possible. For 
instance, if Ni and N,,, are both even and ‘‘m”’ is also even, this can be 
achieved by selecting N,(s) to be the lowest transmission zero below 
the passband and N,,(s) to be the highest transmission zero above. 
The reason for the above requirement is simply the fact that in the 
passband 

N(s)N(—s) 


D(s)D(—8) lente (64) 


will usually be close to a positive constant, and therefore the product 
of (63) and (64) will also have this property, leading to an F'(s)F(—s) 
that is ‘‘small’’ inside the filter passband. 

Finally, we come to the question of separating F(s) out of 
F(s)F(—s). Since F(s) need not be a Hurwitz polynomial and since 
only a few, if any, of the zeros of F(s)F(—s) will be purely imaginary 
(and of even multiplicity), we will have a finite number of possibilities 
to choose from. At this time we have no guideline to offer, we simply 
note that in many cases we had to select roots with alternating real 
parts in order to reduce the size of the coefficients of F(s), as this pro- 
cedure seemed to be necessary to insure the positiveness of the final 
coefficient values. 

Clearly, the selection of F(s) will affect the final network sensitivity 
properties, but further study is needed to clarify the role F(s) plays 
in influencing the sensitivity. 
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APPENDIX C 
Modifying the Feedback Coefficients 


As mentioned in Section I of this paper, the restriction of all feed- 
back coefficients to —1 does not restrict the generality of our structure 
and can be modified to scale the individual 7;(s) blocks. Such a 
scaling is necessary to obtain the maximum overall dynamic range for 
the structure and can be performed simply as follows. 

The 7, transfer function can be multiplied by a constant a, if all 
subsequent transfer functions are modified as 


T r42i 2 an n42i 


#=1,2,---. 
T e42i-1 =a T 421-1 
Qk 


Depending on the number of blocks following 7,, this will either 
leave the overall transfer function unchanged, or will multiply it by 
a,. Also, the —1 in the (kK — 1)th feedback loop must be replaced by 
—1/a,. This procedure works for all blocks. In the case of the first 
block, an alternative is to multiply 7; by a constant ai, multiply the 
second feedback loop coefficient and divide 7’, by the same a: constant. 
This way one can scale each and every one of the 7; transfer functions 
without any effect on the overall transfer function save for a constant 
multiplier. 

For the purpose of adjusting the dynamic range of the individual 
blocks, one must be able to calculate the input and output voltage 
levels of each block. In the case of all —1 feedback coefficients, this 
can be easily done by recognizing that these voltages are numerically 
equal to the branch “voltages” and “currents” of our pseudo-ladders 
of Figs. 6a or b. As such, these can be readily computed by the use of 
continuants again or by any other convenient way. 


APPENDIX D 
Comments on the Positiveness of the Coefficients 


Consider now the realizability of the multiple-feedback structure. 
The synthesis technique given above guarantees the existence of a set 
of T, biquadratic transfer functions with the specified properties for 
any 0< ? S k2,, and arbitrarily specified zero sequence. However, 
some of the coefficients may turn out to be negative. In such a case, 
it is preferable that one or more of the 7; blocks have an overall 
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Fig. 10—Passive bandpass ladder filter. 


negative sign rather than blocks having coefficients of mixed signs. 
We have so far been able to achieve this in all of the many examples 
calculated, but we had to put up with negative 7, blocks in certain 
cases. In particular, nonbandpass filters are likely to have negative 
blocks in their realization. 

However, this is not a shortcoming of the method, since it appears 
even if the structure was derived from a passive double-terminated 
ladder with all positive elements. The sensitivity improvement is still 
realized and no instability will be generated, since our structure is 
minimal. Consequently, there are no pole-zero cancellations and, since 
the overall system poles are the zeros of D(s), the system will be 
stable if the original requirements called for a stable transfer function. 

As an example, consider the passive bandpass filter of Fig. 10. To be 
specific, let us select a passband from 660 Hz to 980 Hz with 0.25-dB 
loss ripple and transmission zeros at 


Zi. = 0 
23,4 = © 
45,6 = + 71280 Hz 
27,8 = ©. 


The resulting poles normalized to the geometric center frequency are 


TABLE [TV—ELEMENT VALUES OF CIRCUIT IN Fia. 10 


Ri = 1.0 C2 = 5.46300 
L; = 0.256378 a 

Gl = 3:71840 Ly = 0.173352 
Lz = 0.981823 Cs = 72.1868 

Co = 0.402048 Le = 0.0129235 
Ch = 0.556146 Re = 0.00327825 
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TABLE V—COEFFICIENTS OF A MULTIPLE-FEEDBACK REALIZATION 
OF THE PASSIVE FILTER SHOWN IN Fia. 10 


1 2 3 4 
No — 1.0 —1.0 0.413433 
i 0.256378 0.394775 — — 
do 1.0 1.98114 5.87644 1.0 
dy 0.256378 — — 0.236646 
de 0.953318 1.76393 5.17359 0.932908 
located at 

Py. = — 0.041320302 + j0.81446903 

P34 = — 0.096338325 + j0.94746237 

Ps,.5 = — 0.085145570 + 71.1203830 

Pz,3 = — 0.030860810 + 71.2239735. 


This filter can be realized as shown in Fig. 10 with element values 
(normalized to 1-ohm input impedance level and 1 rad/s center fre- 
quency) given in Table IV. 
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Fig. 11—Sensitivity improvement. 
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One can derive an active multiple-feedback equivalent of the form 
of Fig. 3 for this structure by writing down the branch equations as 
we have done in Section I, and then eliminate the variables V3 and Vs 
from them. 

The resulting (normalized) coefficients are given in Table V and the 
most noteworthy result is that the third block has a negative sign as- 
sociated with it. On the basis of Appendix C, this is equivalent to 
having a +1 (positive) feedback coefficient in the second and third 
loops. 

Our synthesis technique was then used to generate other (equiva- 
lent) realizations, but we have failed to find one with all positive 
coefficients. Nevertheless, the realization shown above is quite satis- 
factory; the sensitivity improvement compared to the cascade reali- 
zation is shown in Fig. 11. 

Under what conditions are we to accept one or more negative blocks 
in the realization is a question that remains to be answered. 
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Statistical Behavior of Rain Attenuation* 


By S. H. LIN 
(Manuscript received October 25, 1972) 


Thirty-one sets of experimental data on the statistics of microwave rain 
attenuation at frequencies above 10 GHz, in the U.S.A., England, Japan, 
Italy, and Canada, indicate that: (2) the distribution of rain attenuation a, 
in dB, is approximately lognormal with a standard deviation oa of logioe 
ranging from 0.46 to 0.71 for earth-space paths, and from 0.88 to 0.86 
for terrestrial paths; (ii) the distribution of the rain fade duration + ts also 
approximately lognormal with a standard deviation a, of logio 7 ranging 
from 0.44 to 0.76 for both earth-space paths and terrestrial paths. We 
propose a theory to explain this general behavior. A theoretical upper 
bound for the fade duration distribution in the tail region is also given. 

The findings in this paper simplify the determination of rain attenua- 
tion statistics needed for the design of earth-satellite radio links and ter- 
restrial radio links. 


I. INTRODUCTION 


The statistics of rain attenuation are important for the design of 
both terrestrial and earth-satellite radio links using frequencies above 
10 GHz. Many experiments!-* have been performed to obtain data 
on rain attenuation for different frequencies, path lengths, and geo- 
graphical locations. This paper presents the behavior of rain attenua- 
tion statistics found in our study to be common to the available ex- 
perimental data.! 


II. DEFINITIONS 


Let 
V(t) be the time-varying amplitude of the received signal 
voltage normalized to its nonfaded level, 
* An excerpt of this paper has been presented at the 1972 IEEE International 
Conference on Communications at Philadelphia, and included in the Proceedings of 
the Conference. 


t These data are actual measured microwave rain attenuations and not the ap- 
proximate attenuations calculated from rain rate data. 
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a(t) = — 20 logio V (é) be the time-varying rain attenuation 
in dB, 
P(a 2 A) be the expected fraction of time that a(t) exceeds any 
specified value A, 
P., be the expected fraction of time that rain falls at the 
location of the radio link, 
P.(a 2 A) be the expected fraction of raining time that a ex- 
ceeds A, 
Qm be the median value of a during the raining time, i.e., 
P.(a 2 am) = 0.5, 
7(A) be the duration of rain attenuation fades with a ex- 
ceeding any specified threshold A in dB, 
P[r(A) 2 db] be the probability that the fade duration 7(A) ex- 
ceeds any specified duration b, and 
7(A) be the average duration. 


Notice that 
P(a 2 A) = P,:P.(a 2 A). (1) 


Thus, P. may also be called the probability of rainfall and P.(a 2 A) 
the conditional distribution of a under the condition that the rain is 
falling. 


III. SUMMARY OF RESULTS 
3.1 Attenuation Distribution 


(t) The available experimental data on both earth-space paths 
and terrestrial paths in the U.S.A., England, Japan, Italy, and 
Canada show consistently that the conditional distribution, 
P.(a 2 A), is approximately lognormal within the attenua- 
tion range, 1 dB S a S 50 dB, of practical interest.* 

(tt) We propose a theory to explain the lognormal behavior of 
attenuation a(t). In essence, the value of the rain attenuation, 
a(t), at any time instant can be multiplicatively affected by a 
large number of random time-varying parameters of the 
environment such as the present states and the past histories 
of the weather conditions at various locations all over the 
world. The large number of random multiplicative components 
and the central limit theorem lead to the lognormal distribu- 
tion of a(t). 


* This means the distribution of signal amplitude V(t) is approximately 
log-lognormal. 


RAIN ATTENUATION 559 


(171) These results indicate that three parameters, P., o. (the 
standard deviation), and am (the median), are sufficient to 
determine the rain attenuation distribution P(a 2 A). These 
parameters depend on geographic locations. 

(itv) cm increases almost linearly with the path length because the 
median rain rate usually is small and is almost uniform over the 
entire path. 

(v) The dependence of am on frequency follows the theoretical 
prediction of Sezter?° for small rain rate because the median 
rain rate is usually small. The rain rate data in Refs. 13 and 
34 to 37 indicate that the median rain rate ranges from 0.5 to 
10 mm/h, depending on location. 

(vt) oa decreases slightly as the path length increases because of 
the averaging effect of the propagation volume.?*:5 

(vit) oq decreases slightly as frequency increases. 


3.2 Fade Duration Distribution 


(t) The experimental data indicate that the fade duration dis- 
tributions, P[7(A) 2 0], are also approximately lognormal. 
The physical reason for the lognormal duration is the same 
as that for the lognormal attenuation. 

(27) The probability of occurrence of long fade durations has an 
absolute upper bound given by 


1 l X\? 
P (: > x) < -erfe ( eeu ) (2)\* 
2 M 


FZ 





for any X = 1. For example, P(7/7 2 10) S$ 0.0161 means 
that no more than 1.61 percent of the total number of fades 
will have durations longer than 10-7. 


IV. EXPERIMENTAL DATA 
4.1 Rain Attenuation Distribution 


Thirty-one sets of experimental data on rain attenuation distribu- 
tion, each with a time base of six months or longer, are summarized 
in Tables J, II, and III for earth-space paths, long terrestrial paths 
(> 10 km), and short terrestrial paths (S$ 10 km), respectively.' 


* The constant M is defined in eq. (13). 

t In the literature, there are many other sets of experimental data with a time 
base less than six months. Those data are not included because the short-term dis- 
tributions of rain attenuation are fairly random. 


TaB.LE I—Rain ATTENUATION EXPERIMENTS ON EARTH-SPACE PATHS 


Reference Frequency | Path Length Om 
Authors Number Location (GHz) (km) P, oa | (dB) Time Base 

Wilson and Ruscio | 4, 5,17 | Crawford Hill, N. J. 30 Earth-Space | 5.23 X 1077 | 0.57] 1.6 Dec. 8, 1967— 
Feb. 28, 1969 

Wilson and Ruscio | 4, 5,17 | Crawford Hill, N. J. 16 Earth-Space | 5.23 X 107? | 0.71] 0.33 Dec. 8, 1967— 
Feb. 28, 1969 

Wilson 15 Crawford Hill, N. J. 16 Earth-Space | 5.9 XX 107? | 0.67| 0.3 April 1, 1969- 
Aug. 7, 1969 

Wilson 15 Sayreville, N. J. 16 Earth-Space | 5.9 X 107? | 0.60} 0.48 April 1, 1969— 
Aug. 7, 1969 

Wilson 15 Parkway, N. J. 16 Earth-Space | 5.9 XX 107? | 0.70} 0.38 April 1, 1969- 
Aug. 7, 1969 

Davies 31 England 19 Earth-Space | 8.5 XX 107? | 0.46} 0.4 June 1968— 
May 1970 


Note: This table does not include the results of an earth-space path at Point Reyes (30 miles north of San Francisco), California 
by K. O’Brien because the rain attenuation exceeded 8 dB only once (for 8 minutes) during the 9-month period of the experiment at 
that location. 
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TaBLeE IJ—Rain ATTENUATION EXPERIMENTS ON LONG TERRESTRIAL PATHS 


Authors 


Hathaway and Evans 


Zimmerman 
Zimmerman 
Funakawa, et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Turner, Easterbrock, 
et al. 


Blevis, et al. 
Stracca 


Stracca 





Refer- 
ence 


Number 


Location 
Mobile, Ala. 
Mobile, Ala. 


Southern England 


‘Southern England — 
‘Southern England — 
‘Southern England — 
‘Southern England — 


Southern England 








Fre- Path 
quency | Length Am 
(GHz) | (km) P. oa | (dB) 
11 43.5 | 3.7 X 10°? | 0.40; 7.1 
17 22.9 | 3.7 X 107 | 0.46} 6.4 
17 14 3.7 X 107? | 0.49} 4 
12.62 80 10 X 107 | 0.42} 3.1 
ll 24 8.5 X 10° | 0.51 | 0.21 
18 24 8.5 X 1077 | 0.51 | 0.55 
36 24 8.5 & 107? | 0.52] 2.3 
11 24 6.4 X 10°? | 0.40] 0.76 
il 30.6 | 6.4 X 107 | 0.34] 1.2 
ll 38.62 | 6.4 * 10 | 0.38; 1.0 
ll 58 6.4 X 107? | 0.87] 1.4 
15 15.78 | 5.8 X 10° | 0.67 | 0.4 
18 20 6 X10? /0.56) 1.82 
11 20 6 x 10?/0.71} 0.45 








Time Base 
1956 


Aug. 1, 1958- 
Aug. 31, 1959 


Aug. 1, 1958— 
Aug. 31, 1959 


Nov. 1964-— 
Oct. 1965 


Jan. 1967— 
Dec. 1968 


Jan. 1967— 
Dec. 1968 


Jan. 1967— 

Dec. 1968 
1964 

1963 and 1964 
1964 


1964 


May-Oct. 1965 
May-—Oct. 1967 
May-—Oct. 1967 


NOILVONGLLVY NIVU 


T9S 


TABLE [J]—Rain ATTENUATION EXPERIMENTS ON SHORT TERRESTRIAL PATHS 


Authors 


Barnett and Bergmann 
Barnett and Bergmann 


Barnett and Bergmann 


Semplak 
Semplak 
Semplak 
Semplak 


Delange and Dietrich 


Gray 
Kenny 
Hickin 


Refer- 
ence 
Num- 
ber 
































Fre- 
quency 
(GHz) 





Path 
Length 
(km) 


Location 


Palmetto, Georgia 


Palmetto, Georgia 


Crawford Hill, N. J. 


Crawford Hill, N. J. 
Crawford Hill, N. J. 
Crawford Hill, N. J. 


Crawford Hill, N. J. 
Crawford Hill, N. J. 
Merrimack Valley, Mass. 


England 


Time Base 
Nov. 1970-June 1971 
Nov. 1970-June 1971 
Nov. 1970-June 1971 
Jan. 1—Sept. 25, 1970 
Jan. 1-Sept. 25, 1970 
1967 + 1968 + 1969 

1968 + 1969 
1970 
1970 
April 1970-July 1971 


March 1964—Feb. 1966 


Oa 





0.57 
0.6 


0.57 
0.86 

















0.78 
0.86 
0.85 
0.47 
0.55 
0.44 
0.85 























(4B) P, 

0.55 |5.7 X 107 
05 (5.7 X10 
0.55 18.7 xX 107 
0.074 [4.03 x 107 
0.24 | 4.03 X 102 
0.11 | 5.08 X 10 
0.11 | 5.55 x 107 
0.94 | 4.15 x 107 
0.46 | 4.15 X 107 
tl 14 X 107 
0.045| 6.6 X 107 





G9G 
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60 
(1) 30 GHz, NEW JERSEY 

50 (2) 16 GHz, NEW JERSEY 

40 (3) 19 GHz, ENGLAND 


ZA 


PERCENT OF RAINING TIME THAT a@ 





1 10 100 
RAIN ATTENUATION A IN dB 


Fig. 1—Lognormal distributions of rain attenuation on earth-space paths. 


The original data are all given in terms of the unconditional dis- 
tribution P(a# = A). Based upon P, (either available or estimated), 
we convert these unconditional distributions into the conditional dis- 
tribution P.(a = A) by eq. (1). This conversion is done to exclude the 
dry periods in which the rain attenuation is identically zero.* 

When plotted on a lognormal coordinate system, these conditional 
distributions, P.(a 2 A), are all approximately straight lines within 
the attenuation range, 1 dB Sa S 50 dB, of practical interest. 
Figures 1 to 3 show nine examples. The equation describing the 
lognormal distribution is 


(3) 


1 A — a 
Pilg SA). = Pere [=| 


V2.4 


where erfe (~) denotes the complementary error function; o. is the 
standard deviation of logio a during the raining time; and 


Ma = logio Am (4) 


is the mean value of logioa during the raining time. The estimated 
values of c, anda, are given in Tables I, II, and III. 
However, the accuracy of these estimated values of o. and am is 


* Absorption by the clear atmosphere is not considered in this paper. 
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80 


(1) 11 GHz, 43.5 Km, ALABAMA 
70 (2) 12.62 GHz, 80 Km, JAPAN 
(3) 15 GHz, 15.78 Km, CANADA 
(4) 18 GHz, 24 Km, ENGLAND 


60 
50 
40 
30 


ZA 


20 


PERCENT OF RAINING TIME THAT @ 
rS) 


0.1 





0.01 
1 10 100 
RAIN ATTENUATION A IN dB 


Fig. 2—Lognormal distributions of rain attenuation on long terrestrial paths. 


limited by two problems: 


(i) The time base of the experiment may not be sufficiently long 
to yield stable statistics, and 
(tt) some of the published experimental data do not provide the 
probability P, of rain during the experiment.* 
Let 
a(t) = Llogio a(é) J Ha 5) 
Ca 
If a is lognormally distributed, then 6 will be normally distributed with 
zero mean and unity standard deviation. This allows us to pool all 
the available data of 8 on the same graph paper for comparison. Figures 
4, 5, and 6 show the pooled data of the earth-space paths, long ter- 
restrial paths, and short terrestrial paths, respectively. It is seen that 
these experimental results of 6 are indeed normally distributed with 
zero mean and unity standard deviation. 
*The P, values at Crawford Hill, New Jersey, are provided by D. C. Hogg and 
R. A. Desmond from their rain gauge records. The P, values in Alabama, England, 


Japan, and Canada are estimated from the information in Refs. 6, 8, 11, and 13. 
The P, values in Italy and Massachusetts are assumed values using some judgment. 
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30 (1) 18 GHz, 5 Km, GEORGIA 
(2) 18 GHz, 10 Km, ENGLAND 


2A 


PERCENT OF RAINING TIME THAT @ 





1 10 100 
RAIN ATTENUATION A IN dB 


Fig. 3—Lognormal distributions of rain attenuation on short terrestrial paths. 


However, Fig. 6 shows that the data of short paths have significant 
deviations from the normal distribution in the tail region (68 = 2.8). 
Furthermore, most of the deviations in the tail region are downward 
from the straight line approximation. The reason for this nonsym- 
metric deviation is discussed in Appendix B. 






POOLED DATA OF S!X EXPERIMENTS 
ON EARTH — SPACE PATHS 
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Fig. 4—Normal distribution of 8 of earth-space paths. 
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95 
POOLED DATA OF FOURTEEN EXPERIMENTS ON LONG PATHS 
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B eee onan 


2 = i) 1 2 3 4 
Fig. 5—Normal distribution of 8 of long terrestrial paths. 
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POOLED DATA OF ELEVEN EXPERIMENTS ON SHORT PATHS 
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Fig. 6—Normal distribution of 8 of short terrestrial paths. 
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(1) 2.6 Km, CRAWFORD HILL, N.J. 

(2) EARTH — SPACE, CRAWFORD HILL, N.J. 
3 (3) 20 Km, ITALY 

(4) 24 Km, ENGLAND 


Qm IN dB 


0 10 20 30 40 
FREQUENCY IN GHz 





Fig. 7—Effect of frequency on median attenuation am. 


4.2 Dependence of o, and am on Path Length and Frequency 


The experimental results, discussed in Section 4.1, show that the 
three parameters, P., oa, and am are sufficient to determine the dis- 
tribution P(a 2 A). Therefore, it is important to study the depend- 
ence of o, and a» on path length and frequency. 

The effects of path length and frequency on o, and a» are shown in 
Figs. 7 to 10. These experimental results indicate the effects of path 


1.5 
(1) 2.6 Km, CRAWFORD HILL, N.J. 
(2) EARTH — SPACE, CRAWFORD HILL, N.J. 
(3) 20 Km, ITALY 
(4) 24 Km, ENGLAND 
1.0 


(Ie 





0 10 20 30 40 
FREQUENCY IN GHz 


Fig. 8—Effect of frequency on standard deviation oz. 
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(1) 18.5 GHz, CRAWFORD HILL, N.J. 
(2) 17 GHz, MOBILE, ALABAMA 
(3) 11 GHz, ENGLAND 


Qm IN dB 





0 10 20 30 40 50 60 
PATH LENGTH IN Km 


Fig. 9—Effect of path length on median attenuation am. 


length and frequency on oa and am as stated in terms (2v) to (viz) of 
Section 3.1.* 

The theoretical calculation®® indicates that the increase of rain at- 
tenuation with frequency f is slightly faster than the square law in the 
range: 10 GHz Sf S 60 GHz, and 1 mm/h S rain rate S$ 15 mm/h. 
This is the basis of the dashed curves in Fig. 7. The slopes of straight 
lines in Fig. 9 are proportional to median rain rate and extinction 
coefficient. Curves (1) and (2) in Fig. 9 indicate that the median rain 
rate in Alabama is much larger than that in New Jersey. 


4.3 Fade Duration Distribution 


The available nine sets of experimental data on the histogram of the 
durations 7(A) of rain attenuation fades are summarized in Table IV.t 
We convert these histograms into the cumulative distribution 
P[7r(A) 2 6]. On a lognormal coordinate system, these fade dura- 
tion distributions are all approximately straight lines. Figure 11 shows 
two examples. 

* In Fig. 8, the ¢ in England seems to increase slightly with frequency in contrast 
to those in the U.S.A. and Italy. A possible reason for this inconsistency is that the 
time bases for the three sets of data for 11, 18, and 36 GHz measured in England 
are not concurrent. 


_'Some of the available data in the literature are not included because of short 
time base. 


RAIN ATTENUATION 569 


1.5 


(1) 18.5 GHz, CRAWFORD HILL, N.J. 
(2) 17 GHz, MOBILE, ALABAMA 
(3) 11 GHz, ENGLAND 





ie) 10 20 30 40 50 60 
PATH LENGTH IN Km 


Fig. 10—Effect of path length on standard deviation og. 


The equation describing the lognormal distribution of 7(A) is 


logo b a | 
Vlo; 


where yu, and a, are the mean and the standard deviation respectively 
of logio (A). The estimated values of o, and the average fade dura- 
tion 7(A) are given in Table IV. Again, the unstable statistics, caused 
by insufficient time bases, limit the accuracy of these estimated values 
of o, and 7. 


P[r(A) 2 6] = perf (6) 


80 


(1) 30 GHz, 21 dB FADE THRESHOLD, N.J. 
(2) 19 GHz, 5dB FADE THRESHOLD, ENGLAND 


b 


70 


60 
50 


40 
30 
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20 


PERCENT OF FADES FOR WHICH 7 





1 10 100 
FADE DURATION b IN MINUTES 


Fig. 11—Lognormal distributions of fade duration. 


Authors 
Barnett 


Wilson and Ruscio 
Wilson and Ruscio 


Wilson and Ruscio 


Davies 


Semplak and 
Turin 


Semplak and 


Turin 


Semplak and 
Turin 


Stracca 


TABLE [V—EXPERIMENTS ON DISTRIBUTION OF FADE DURATION 


Refer- 
ence 
Number 


Location 


Frequency 
(GHz) 


Path 
Length 
(km) 


Or 


T(A) 


(min) 


Fade 
Thresh- 


Time Base 


November 1970—June 1971 
Dec. 2, 1967—-Feb. 28, 1969 
Dec. 2, 1967—-Feb. 28, 1969 
Dec. 2, 1967—Feb. 28, 1969 
June 6, 1968-May 31, 1970 
June 23, 1967—Oct. 31, 1967 


June 23, 1967-Oct. 31, 1967 


June 23, 1967—Oct. 31, 1967 


May-Oct 1967 


OLS 


S161 TlddV “IVNUNOL TIVOINHOGDL WALSAS TIAG THL 


RAIN ATTENUATION 571 


Vv. THEORY 
5.1 Rain Attenuation Distribution 


Existing theory!?8 and experimental data!~99-%3 indicate that 
rain attenuation is a complicated function of many parameters of the 
propagation medium: the total number of rain drops in the path, the 
drop size distribution, the fine grain spatial characteristics of the rain 
density along the path, wind velocity, the presence of up or down 
drafts, raindrop shape, raindrop cant angles, the storm cell shapes and 
sizes, raindrop temperature, etc. In other words, the attenuation a(t) 
is a function of many random time-varying parameters of the medium. 
Furthermore, through the coupling of the atmosphere, the above- 
mentioned parameters of the propagation medium depend on the 
present states and the past histories of the weather conditions at 
many near or faraway locations and altitudes. 

We will assume that the rain attenuation can be affected by a large 
number of random time-varying multiplicative components: 


a(t) = Sit) -S2(t)-Ss(t)- + Sr). (7) 


Each of {S;(t)}{=% represents the random modification factor due to 
an environmental parameter. 
Taking logarithms on both sides of eq. (7) yields 


loga = log 8S: + log S2 + --- + log S, (8) 


which shows that log a is a summation of a large number of random 
variables. Then by the central limit theorem,®* the distribution of 
log a approaches a normal distribution for large n if there is no domi- 
nant component. Therefore, the distribution of a is approximately 
lognormal.*t 
The basis for the multiplicative formulation (7) is that the environ- 
mental parameters affect the rain attenuation a(t) in a proportional 
fashion (i.e., in terms of percentage) rather than an additive fashion. 
For example, at 30 GHz frequency, the theoretical calculation shows 
that, when the rain temperature decreases from 20°C to 5°C, the rain 
attenuation increases by approximately 4 percent. This means the 
*It is interesting to note that, in Fig. 12 of Ref. 10, the lognormal distribution 
appears to be a reasonable fit to the probability, distribution of rainfall rate. The 
physical reason for this behavior probably is similar to that for rain attenuation a 
disciiesed in this section. 
A basic characteristic of a lognormal random variable y is that its value can vary 
in the entire semi-infinite range: 0 S$ y <. For the rain attenuation problem, 
the attenuation a(dB) possesses this basic characteristic. On the other hand, the 


normalized signal amplitude V (¢) is strictly confined to the finite range: 0 S V(t) < 1, 
which rules out the possibility of a lognormal distribution for V (é). 
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variation of rain attenuation, due to 15°C variation of rain tempera- 
ture, can be 4 dB, 0.4 dB, or 0.04 dB if the attenuation at 5°C is 100 
dB, 10 dB, or 1 dB, respectively. Similar arguments apply to the 
effects of other parameters on the attenuation. A more general inter- 
pretation of formulation (7) is discussed in Appendix A. 

We emphasize that some of the components {S;(t)}'=? may be 
extremely slowly varying functions, which may take several months 
or even several years in order to show their effects. For example, the 
rainfall intensity-duration-frequency data of the Weather Bureau*:* 
show that in New Jersey the return period” for a rain rate exceeding 
150 mm/h, which continues for a 5-minute duration, is about 5 years. 
Therefore, if the time base of a rain attenuation experiment is less 
than 5 years, the chance of missing these rare and extreme events is 
very high. The justification for considering such a long-term distribu- 
tion of a(t) is that the microwave radio systems, which are designed 
based on these statistics, contains many repeaters, each sampling its 
own rain universe and contributing to total path outage. 


5.2 Fade Duration Distribution 


We also assume that the duration 7(A) of a rain attenuation fade 
with a 2 A is affected by a large number of random time varying 
multiplicative components 


T(A) = X Xo: X3:+ Xn. (9) 


Each of {X,}=7 represents the random modification factor of an 
environmental parameter. Therefore, the long-term distribution of 
7(A) is also approximately lognormal. 


5.3 Upper Bound for Fade Duration Distribution 


In radio system design, one is concerned with the occurrence prob- 
ability of an unusually long continuous outage. It is desirable to have a 
“quick estimate” of the fade duration distribution, especially in the 
tail region of long duration. In our experimental and theoretical 
study“!-® of the lognormal distributions of durations of rain and 
multipath fading, W. T. Barnett*! has found an upper bound on the 
fade duration distribution as discussed in the following. 


Let* 


z= 


aia 


(10) 


* The idea of normalizing fade duration r.to the average duration 7 originates 
from the work’ of S. O. Rice on fade duration distributions. 
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The lognormal] distribution of x can be written as 


logio xX — Le 
eae 
where yw and o are the mean and the standard deviation, respectively, 
of logio x. For lognormally distributed x, it is easily shown that 

E = eb! M+} (02/ M2) (12) 


Pia 2X) = jefe] (11) 


where Z is the mean value of 2, and 
M = logie = 0.484. (13) 
The definition (10) implies that 


T 
#=-=1. (14) 
T 
Equations (12) and (14) show that 
_ 1 o? (15) 
"9M 


Substituting (15) into (11) yields 


logio X + aS} a 


P(x = X) = derf 
(2 2 X) s erfe| Bo 


Therefore, the lognormal distribution of x is completely determined 
by only one parameter co. 


N 
(=) 


_— — — Pmax (= B10)=1.61% 


P( 2210) IN PERCENT OF NUMBER OF FADES 





oO 0.5 1.0 1.5 2.0 2.5 3.0 
STANDARD DEVIATION o& 


Fig. 12—Effect of standard deviation ¢ on the probability of long duration of fade. 
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which implies that for any X = 1 the probability P(t = X) as a 
function of o has a maximum value: 


logio xX 
Piste sos perfe| 7 | (18) 


o = V2M logis X. (19) 


For example, Fig. 12 shows P(x = 10) as a function of o. 


Figure 13 shows the maximum probability Pmax(v 2 X) as a fune- 
tion of X as given by eq. (18). By the use of “‘C-discriminant equa- 
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Fig. 14—The envelope of the family of lognormal duration distributions. 
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tion” in the theory of ordinary differential equation,* it can be shown 
that eq. (18) is, in fact, the envelope of the family of lognormal dis- 
tributions (16) with o as the family parameter. In other words, the 
family (16) are all tangent to (18) as shown in Fig. 14. The shaded 
area in Fig. 13 is a forbidden region where the lognormal duration dis- 
tribution will never penetrate. Therefore, the upper bound for the 
fade duration distribution is given by eq. (2). 


VI. CONCLUSION 


Both the experimental data and theory indicate that both the rain 
attenuation distribution and the fade duration distribution are 
lognormal. The detailed results have already been given in Section III 
(Summary of Results). 
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APPENDIX A 
Generalized Interpretation of Formulation (7) 


From a more general viewpoint, the effect of an environmental 
parameter Z,(t) on a(t) may be more complicated than the simple 
linear proportional relation 


a(t) « Z(t), 7= 1,2, +--+, n. (20) 


For example, a() may depend not only on the present value, but also 
on the past history of the environmental parameter Z;(¢); then the 
relation between a(t) and Z;(t) becomes 

t 


a(t) « r H(t, t’)Z;(t’) dt’ ais 


Pl, 2.284 


where H,(¢, t’) is the impulse response of a(t) if the input Z;(¢) is an 
impulse 6(¢ — ¢’) applied at ¢’. In order to allow for the more general 
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and complicated relation between a(t) and Z;(t), we shall use the 
mathematical operator notation: 


a(t) = GiLZ.(t) ], t= 1, 2, 3, see ee (22) 


where G; is a transformation from Z;(t) into a(t). Then the formulation 
becomes 
a(t) = Gi(Z1) *G2(Z2) +++ Gn (Zn). (23) 
Let 
Sid) = G[Z,0] 7 = 1, 2,3, ---,n. (24) 


Substituting (24) into (23) leads to (7). 

The formulation can be further generalized to include the cases 
where the effects of various environmental parameters on @ are not 
completely separable. Let 


vy = logioa, (25) 
ci€; = logio Si = logio Gi(Z;) 

t= 1, 2,3, ---,n—1, and (26) 

R = logio Sn. (27) 


Substituting (25), (26), and (27) into (8) gives 


y (£1, £2, +++, Ena) = 7 + Ci(—s — &) + Co(& — &) + -- 
+ Ca-1(Ea-1 — En) + (R - B (28) 


where 

F¥ = Cyéy + Coke + +++ + Cp-rén-1 + R.z (29) 
Equation (28) can be interpreted as the first-order Taylor series ex- 
pansion of y(£1, £, ©:+, .-1) With a remainder term,* R — R. The 


main point in eq. (28) is that the effects of the (n — 1) random vari- 
ables {£;}:=7~! on y do not have to be completely separable because 
of the remainder term R — R. If R — R does not dominate the sum in 
the right-hand side of eq. (28), then, by central limit theorem, the 
distribution of y (i.e., logioa) is approximately normal even if the 
effects of various environmental parameters are not completely 
separable. (This argument originates from Refs. 46 and 47.) Therefore, 
the formulation (7) includes very general and complicated relation- 
ships between the environmental parameters and rain attenuation a(t). 


* If higher-order derivatives of y(&, &, +++, &n-1) exist, then R — R represents 
the sum of all the higher-order terms; otherwise, R — R represents the difference 
between y(é1, £2, +++, &n-1) and its first-order Taylor series expansion. 
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APPENDIX B 
Deviations of Experimental Data From Lognormal Distribution 


Aside from the experimental error, there are two major factors which 
contribute to the deviations of experimental data from the lognormal 
distribution. 


B.1 Effect of Finite Number of Components 


Davenport and Root**® have indicated that when the number, n, of 
components is finite, the normal distribution may well give a poor 
approximation to the tazls of the distribution of the sum (8) even 
though the limiting form of the sum distribution is, in fact, normal. 
Therefore, in practice, we believe that the deviation of the experimental 
data from the lognormal distribution may increase toward the tails. 


B.2 Effect of Time Base 


Since some of the components in eq. (7) are extremely slowly 
varying, reducing the time base will reduce the number of contribut- 
ing components since the influences of slow components are approxi- 
mately constant in a short experiment. Therefore, we expect the 
deviation of the experimental data from the lognormal distribution to 
increase as the time base decreases. 

The minimum required time base for the convergence of the experi- 
mental data to the lognormal distribution increases as 


(t) path length decreases, or 
(iz) operating frequency decreases, or 
(zit) attenuation range of interest increases. 


I'urthermore, when the time base of a rain attenuation experiment 
is not long enough, the deviations of the short-term distribution 
P.(a 2 A) from the long-term distribution P(a = A) are usually non- 
symmetric in the deep fade region, i.e., Ps(a = A) is more likely to 
be less than P(a 2 A). This nonsymmetry is caused by: 


(i) The distribution of fade duration 7(A) is lognormal, which is 
nonsymmetric with respect to the average fade duration 
7(A). Typically, about 70 percent of fade durations are 
shorter than the average duration. 

(iz) The probability distribution of the number N(A, T) of ob- 
served deep fades, exceeding the margin A dB in the period T, 
is somewhat similar to a Poisson distribution which is also 
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nonsymmetric with respect to the average number N (A, T) of 
deep fades in a period T. 


Since, when the time base is too short, both N(A, 7) and r(A) have 
a higher chance of being less than their average values N (A, T) and 
7(A), respectively, then the short-term distribution P,(a 2 A) tends 
to deviate downward from the long-term distribution in the tail 
region.* 


B.3 Effect of Time Base on Standard Deviation oa 


Since it takes a long time base to include appreciable effects of 
slow components, we expect that co, increases slightly as the time base 
increases, and reaches an asymptotic value only after the time base 
is long enough to include the effects of all the possible slow components. 
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Simple, Low-Loss Joints Between 


Single-Mode Optical Fibers 


By C. G. SOMEDA* 
(Manuscript received September 21, 1972) 


Low-loss joints between single-mode optical fibers have been made with- 
out microscopic alignment, without fusing the tips, and without monitoring 
the transmitted power while the joints are assembled. The fibers are tightly 
held in an embossed groove; an index-matching liquid is added. Average 
power coupling efficiencies close to 90 percent in the red and to 85 percent 
in the infrared have been obtained. Mediocre end faces are acceptable. 
Realistic discrepancies between the fiber cladding diameters (slightly in 
excess of twice the core diameter) do not deteriorate the results. 


I. INTRODUCTION 


Recently, several authors have dealt with techniques for making 
low-loss joints between multimode! or single-mode’ optical fibers. 
For long-distance communication channels, joints between fiber 
cables, simply assembled in the field, would be very valuable. 

The joints that were made so far with the technique described in this 
paper connected just one pair of fibers at a time. However, this is in- 
tended to be the first step in the development of joints between multi- 
fiber cables. Therefore, in the present experiments emphasis has been 
put on getting repetitive results which are not a consequence of careful 
laboratory adjustment. This is an important difference with respect to 
the methods that gave the best results reported so far.!:> Nonetheless, 
the joints between single-mode fibers described here compare favor- 
ably with previous ones as far as the losses are concerned. The demon- 
stration with single-mode fibers (core diameters of ~4 um) means that 
the technique applies also to multimode fibers, the larger core diam- 
eters of which should make the alignment less critical. 


* This work was performed at Bell Telephone Laboratories, Incorporated, Craw- 
ford Hill Laboratory, Holmdel, New Jersey, during a six-month internship sponsored 
by a NATO-CNR fellowship. 
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Fig. 1—View of the unassembled parts of a sandwich joint. 


A complete description of the technique follows in Section II. As 
a rough sketch, we may say that the tips of the fibers are aligned in 
a groove (Fig. 1), embossed in a Plexiglas sheet.* The fiber tips are 
pushed against each other. An index-matching liquid is added. A 
“sandwich” is formed with a flat Plexiglas sample, and then squeezed by 
means of a small vise (Fig. 2). No fusion takes place and no microscopic 
alignment is needed. 

It is obvious that the fiber claddings are aligned in this way and 
therefore the results depend on the fiber core being coaxial with the 
cladding. The experiments reported here were performed with single- 
mode fiberst having a core diameter of 3.7 um and a cladding diameter 
of 254 um; coincidence of the core and cladding axes within ~1 wm 
had been observed preliminarily. 

Experiments were performed at 6328 A and at ~9000 A. Repetitive 
measurements gave average power couplings around 87 percent in the 
red and around 88 percent in the infrared ; best results were 93 percent 
and 87 percent, respectively. These figures were obtained without pay- 
ing any particular attention to the quality of the end surfaces of the 
two fibers. 

The experiments were completed by testing the effects of a difference, 
d, between the cladding diameters, as described in Section IV. It turned 
out that sandwiches of good mechanical quality exhibit no deteriora- 
tion of the results for d= 10 um (i.e., more than twice the core diam- 
eter). For d= 30 um, the results were deteriorated, but still showed 


* Plexiglas (methyl methacrylate), registered trademark of Rohm and Haas. 
These fibers were manufactured by Corning Glass Works, Corning, N. Y. 
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Fig. 2—View of an assembled joint. The length of the visible side of the sandwich, 
perpendicular to the fiber, is ~2 cm. 


the presence of an alignment mechanism, which is believed to be due 
to surface tension in the index-matching liquid. 


II. DESCRIPTION OF THE JOINING TECHNIQUE 


In order to obtain a joint of the type shown in Fig. 1, the first opera- 
tion to perform is to emboss a Plexiglas sheet with a groove that fits 
the dimensions of the fibers to be joined. 

A simple and economical technique for embossing grooves in a 
thermoplastic substrate by means of a glass fiber has been reported by 
Ulrich e¢ al.° A further simplification of this procedure proved success- 
ful; there is no need to heat the substrate. A rectangular Plexiglas 
sample and a fiber piece, the length of which slightly exceeds the Plexi- 
glas size, are put between two milled aluminum blocks (flat within +1 
mil) and then tightly pressed by means of a vise for a few minutes. 
When the vise is released, the fiber separates from the Plexiglas, leav- 
ing a sharply embossed groove. Any other piece of a nominally equal 
fiber can be introduced into the groove and fits it very closely. Figure 3 
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Fig. 3—Microscopic picture of a single-mode fiber tip (on the left) lying in an 
embossed alignment groove. Magnification 200X (1 cm in the picture is 50 um in 
actual scale). 


is a microscopic picture (magnification 200) of a fiber tip lying in 
a groove, embossed with another fiber. Figure 4 is a further enlarged 
view (magnification 500) of the groove edge in the vicinity of the 
fiber end; it shows no defect comparable to the size of the fiber cores to 
be aligned. It also shows a fairly regular pattern of longitudinal 
wrinkles due to the compression of the Plexiglas som 

First, Plexiglas having a nominal thickness of jg inch was embossed 
with quartz fibers, 254 um in cladding aiarieten The samples were 
permanently curved with a remarkable convexity of the embossed side. 
The cover (Fig. 1) then tended to flip. Using a thicker Plexiglas sample 
(e.g., = § inch), the curvature was negligible. 

There is some evidence that the deformation of the embossed Plexi- 
glas sheets is elastic (at least partially) with a long time constant. 
Therefore, the joints have to be made with newly embossed samples, 
or with samples where a fiber has been constantly pressed in. 

To assemble a joint, the fibers are put in the groove, a small amount 
of a suitable index-matching oil is added, and the cover is placed on 
top. At first, the top is not pressed against the bottom. Next, one 
fiber is pushed against the other in the axial direction; finally, the 
sandwich is squeezed with a small vise. 
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Fig. 4—Microscopic picture of the edge of a groove, with a fiber tip on the left. 
Magnification 500X (1 cm in the picture is 20 um in actual scale). 


It seems that the index-matching oil is rather slow in wetting the 
fiber tips, and some time is required for air bubbles to escape. The oil 
also has the purpose of lubricating the fiber longitudinal motion. Some- 
times an excess of it deteriorated the joint performances, probably be- 
cause an unwanted axial movement of the fiber tips was produced by 
the oil flux when the sandwich was squeezed. However, the presence of 
a small amount of liquid is needed in order to get good coupling effi- 
ciency (see Section 3.2). Furthermore, the presence of the fluid is be- 
lieved to be responsible for the excellent behavior in case of unequal 
diameters of the two fibers (Section IV). 

The amount of light transmitted by the joint does not always in- 
crease with increasing vise pressure. Often, though, a decrease in the 
coupling beyond a certain pressure was interpreted as due to some mis- 
alignment caused by careless movements accompanying the vise 
tightening. When the groove was deep enough and the top could not 
swing or slide, then such a decrease in transmitted power was absent 
or negligible. A good rule is to press the joint between surfaces that are 
not too stiff, so that the force is distributed on the whole area of the 
sandwich. 
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III. MEASUREMENTS 


3.1 Tests at 6328 A 


Sandwich joints of the type described in Section II were tested first 
at the wavelength \ = 6328 A. Based on manufacturer’s data, the cor- 
responding normalized frequency,’ V = 2.2, means that single-mode 
propagation takes place. Direct observation of the far field at the end 
of a fiber confirmed that. 

The set-up is illustrated in Fig. 5. The light emitted by a commercial 
He-Ne laser was chopped and then launched into the input fiber by 
means of a microscope objective (magnification 40). A curved part 
of the fiber was lying in an index-matching liquid, to strip off the light 
launched into the cladding. The length of this mode-stripping section 
was such that any further addition or small subtraction did not affect 
the power level at the output of the fiber. This output was detected 
by means of an Si solar cell and the signal was sent to a lock-in ampli- 
fier. Figure 5 shows both arrangements that are needed in order to 
evaluate the insertion loss of a joint: a “‘reference’”’ arrangement, where 
the detector is directly connected to the input fiber, and a ‘‘measure- 
ment”’ arrangement, where a joint and an output fiber are added. 
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Fig. 5—Block diagram of the set-up used for the measurements at 6328 A. The 
broken lines (Part 1) refer to calibration (measurement of a reference level). The 
solid lines (Part 2) refer to the actual measurement of the light transmitted through 
a joint. 
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TABLE I—PoweER CovupLiING EFFICIENCIES MEASURED AT 6328 A ann 
CORRESPONDING INSERTION LOSSES OF THE JOINTS 


























Sample Number 1 2 3 
Best results 92% 93% 92% 
(0.35 dB) (0.3 dB) (0.35 dB) 
Worst results 74% 80% 81% 
(1.3 dB) (1.0 dB) (0.9 dB) 
Average results 80.7% 88.7% 86.6% 
(0.9 dB) (0.5 dB) (0.6 dB) 
Number of measurements 12 6 6 
Number of different ends 4 2 2 











Attenuations measured in this way result from both the joint inser- 
tion loss and the losses in the output fiber. Previous data® and direct 
observation show that the output fiber attenuation is negligible with 
respect to the loss in the joint, despite the need for a length of the out- 
put fiber that could insure an effective stripping of all the light trans- 
ferred by the joint into the cladding. 

This way to perform the measurements was preferred to that where, 
after taking a reference, the fiber is cut and then joined again, because 
frequent checks of the reference were needed. Indeed, launching into 
a single-mode fiber is so critical (the power drops by 3 dB if the fiber 
tip is misaligned by a few microns), that accidental causes can pro- 
duce strong shifts in the reference level. For example, in our case the 
material used to bind the fiber tip on an x-y-z manipulator could deform 
slowly, under the strength applied by the curved and quite elastic 
fiber. This effect was compensated for by optimizing the alignment of 
the fiber with respect to the incoming beam before any meter reading, 
both in the reference and in the measurement arrangements. Discrep- 
ancies between optimized reference levels monitored before and after 
a measurement were never larger than 3 percent and usually much 
smaller than that. When the reference levels monitored before and 
after a measurement were different, the more pessimistic estimate of 
the joint losses was taken. 

Table I contains the best, worst, and average results of three se- 
quences of measurements. As for the ‘worst’? results, they exclude 
only those instances where clearly identified man-made mistakes took 
place; i.e., axially separated fibers, one fiber lying out of the groove, and 
Plexiglas debris shaved by the fibers and accumulated between the two 
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ends. These errors were very rare, and were easily detected because of 
a large amount of scattered light shining from the sandwich. This point 
will be discussed further in Section 3.2. 

The sequences of Table I refer to assembling, several times, joints 
that make use of three embossed grooves, one in a jg-inch-thick plate 
and two in $-inch-thick plates. All the measurements were independent, 
in the sense that after any of them the joint was at least disassembled 
completely and then reassembled. The ends of the fibers were changed 
often, as shown in Table I. The amount of index-matching oil was 
changed quite often too. 

An important point is that none of the measurements were of a joint 
where the two coupled end faces resulted from one cut of a fiber; they 
always resulted from two independent breaks, in order to simulate real 
situations. 

A remarkable advantage of this technique is that there is no need for 
very accurate flatness of the fiber end faces. As long as there are no lips, 
which would prevent the two cores from getting close to one another, 
results are good. Figure 6 is a microscopic picture of one typical pair 
of ends used in the measurements. Attempts to get end surfaces of this 





Fig. 6—Microscopic picture of a typical pair of fiber ends used in the sequence of 
experiments summarized in Table I, The bright zone in the center is not the fiber 
core, but Just a consequence of the illumination. 
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Fig. 7—Effects of the longitudinal separation between the fibers on the coupled 
power. Lines A and B, from Ref. 2, represent measured values (without and with 
index-matching oil, respectively) at 6328 A when the fibers were aligned by means of 
micromanipulators. The vertical segments represent the range of results obtained 
when the fibers were held in a sandwich joint. 


quality were almost always successful if the fibers were pulled after 
scoring them with a lathe tool or with an electrical discharge generated 
by a Tesla coil. Many successful attempts simply consisted of clamping 
the fibers in a Plexiglas sandwich and then pulling them. Consequences 
for the multifiber case will be stated in Section V. 

Another interesting observation was that when a joint was tightened, 
bringing the fiber cores into alignment, the forward scattered light on 
the side surface of the output fiber decreased very remarkably. 

Some additional measurements were performed in order to establish 
the effects of longitudinal separation between the fibers. The results 
are compared with those obtained by Bisbee? in Fig. 7. The very large 
spread of values observed for a 5-mil separation can be justified by an 
occasional shortage of index-matching oil, which could cause an air gap 
sometimes to show up when the fibers were separated ; or by deviations 
of the groove from a rectilinear shape, which would cause unwanted 
transverse displacements of the fiber tips to accompany the longitudi- 
nal ones. The best results (7 = 56 percent for a 2-mil separation, 
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TasiLe [JI—Power Courtine Erriciencies MEASURED AT 9000 A AND 
CORRESPONDING INSERTION LOSSES OF THE JOINTS 








Sample Number 2 3 
Best results 87.5% 87% 
(0.55 dB) (0.6 dB) 
Worst results 78.5% 80% 
(1.05 dB) (1.0 dB) 
Average results 82.5% 83.7% 
(0.85 dB) (0.75 dB) 
Number of measurements 11 13 


Number of different ends 2 2 


n & 27 percent for a 5-mil separation), which are in very good agree- 
ment with Bisbee’s,? were observed more frequently than the bad ones. 

In multifiber joints, it is unlikely that the axial separation of each 
pair of fibers will be as small as it was for the single pairs tested so far. 
However, one can think of enlarging the waist of the light beam in the 
fibers, by using a smaller normalized frequency’; the confocal length of 
the beam would then grow as the square of the beam waist. The effect 
of an axial separation upon the coupling efficiency would be smaller 
than that shown in Fig. 7. 


3.2 Tests at 9000 A 


It was pointed out before that, while working with visible radiation, 
a misalignment or a gap between the fibers was revealed by a large 
amount of scattered light leaking from the joint. In a multifiber joint, 
this light will not identify the pair of fibers that form a leaky connec- 
tion. Hence, it was necessary to check whether the very small number 
of mistakes in assembling the joints and the high average of their per- 
formances were independent of the information that the operator was 
provided by the scattered light. 

A sequence of tests in the near infrared, at \ 9000 A (where 
V = 1.55), was then performed using a small-area GaAs LED.* The 
set-up differs from Fig. 5 only because the chopping of the signal is per- 
formed by an audio oscillator driving the LED. Launching spatially 
incoherent radiation into the fiber is much less critical than in the case 
of the laser; this results in a much smaller long-term drift of the refer- 
ence level. 

* The LED was built and provided by C. A. Burrus.® 
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Tests were performed on joints that made use of the grooved samples 
referred to in columns 2 and 8 of Table I. Quite repetitive results were 
obtained. They are summarized in Table II. Comparison with Table 
I shows a very slight deterioration of the best results (power coupling 
coefficients from 92-93 percent down to 87-87.5 percent) but even 
smaller changes of the averages (from 80-81 percent down to 78.5-80 
percent) and almost no change in the lower ends of the range. Let 
us re-emphasize that the tested joints were not assembled while moni- 
toring the transmitted power. 

Comparison of both Tables I-and II with previously published data 
on single-mode fiber joints® shows the desirability of the present sand- 
wich technique. Best results given by the fusion technique® show power 
transmission efficiencies of ~80 percent, at a normalized frequency 
V = 2.0; the average results given by the sandwich technique are 
better, both above and below that value of V. 

The number of times when the attempts to assemble the joints were 
unsuccessful because of easily identified mistakes in the fiber position 
increased to three in a sequence of twenty-seven measurements, com- 
pared to one in twenty-five measurements at 6328 A. This value is still 
low enough to insure that the technique is suitable for further develop- 
ments in a multifiber system, where the risk of making such mistakes 
has to be negligible. 

Some information on the importance of the index-matching oil has 
been collected, too: “dry” joints gave best power coupling efficiencies 
in the order of 50 percent. Comparison with available data? shows that 
the oil produces some other effect besides the elimination of end reflec- 
tions. Quite likely, it compensates for irregularities of the end faces; 
also, it provides some alignment mechanism via surface tension. The 
last point is discussed further in the next section. 


IV. EFFECTS OF DISCREPANCIES BETWEEN THE FIBER OUTER DIAMETERS 


The usefulness of the sandwich technique would be limited if its 
performances were sensitive to discrepancies between the diameters of 
the two joined fibers, a situation that will occur in practice. The results 
reported in this section show that the sensitivity to this kind of imper- 
fection is small enough to satisfy practical purposes. 

A controlled difference, d, between the outer diameters of the fibers 
was introduced by etching one fiber with commercial hydrofluoric acid. 
First, a value d = 10 um was chosen, because it would be sufficient to 
allow more than a complete offset of the two cores if the claddings were 
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TasBLteE III—Power Covurpiine EFFICIENCIES OF JOINTS BETWEEN 
FIBERS WITH A 10-um DIFFERENCE IN CLADDING DIAMETERS. 
MEASUREMENTS PERFORMED AT 9000 A. 








Best result 86% (0.65 dB) 
Worst result 82% (0.85 dB) 
Average result 84.7% (0.7 dB) 
Number of measurements 7 

Number of different ends 2 





aligned along one generatrix; besides which, 10 um is 4 percent of the 
cladding diameter, which can be thought of as a realistic tolerance. 

A new sequence of tests was performed at \ = 9000 A, using the 
embossed sample referring to column 2 in Tables I and II. The results 
are summarized in Table III. Comparison with Table II shows that 
there is no appreciable deterioration of the joint performances. On the 
other hand, it was observed that the transmission was sometimes (not 
always) sensitive to small changes in the force applied by the vise; 
e.g., for one of the best joints assembled without monitoring the trans- 
mitted power (7 & 86 percent), small changes in the tightening could 
cause a drop of the transmitted power to 82 percent or a raise to 89 
percent (insertion loss ~0.5 dB). 

In all cases, releasing the vise slowly caused, first, a decrease in 
transmitted power; then a maximum, close to the figures quoted be- 
fore, showed up; later, the signal decreased monotonically down to the 
noise level. The secondary maximum, which was very sensitive to vibra- 
tions and shocks, is believed to be entirely due to an alignment caused 
by the surface tension of the index-matching oil. The same phenomenon 
is believed to contribute to the excellent alignment under the strongly 
tightened vise, leading to the figures reported before. 

Later, the difference between the cladding diameters was brought to 
a value d& 30 um, i.e., about 12 percent of the outer diameter. The 
following behavior was observed: 


(1) There was practically no coupling as long as the joints were not 
tightened very strongly. This was a remarkable difference with 
respect to the case of equal diameters, where the simple contact 
of the two faces, before squeezing the sandwich, could often 
produce a power transmission in the range 10 to 30 percent. 
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(it) Several joints, assembled and very tightly squeezed without 
monitoring the transmitted power, had coupling efficiencies up 
to 40 to 50 percent and were not very sensitive to vibrations 
and shocks; however, the worst results ranged down to ~10 
percent (very seldom) or 20 percent (more frequently). 

(117) Repeated operations of tough tightening followed by partial 
releasing, while the transmitted power was monitored by a 
meter, led to coupling efficiencies in the range 80 to 84 percent. 
These joints were not tightly squeezed, and were very sensitive 
to vibrations and shocks, which could reduce the coupling all 
the way to 10 percent. Usually, though, the joints that had 
been disturbed, once abandoned to themselves, tended to re- 
store spontaneously power coupling efficiencies up to 70 to 75 
percent. 


These observations and comparison with Bisbee’s data? show that 
a remarkable alignment mechanism, due to the presence of the liquid, 
is still active for such a large difference between the cladding diameters, 
even if the performances are not any more suitable for practical 
applications. 


Vv. CONCLUSION 


Sandwich joints between single-mode fibers, aligned in embossed 
grooves, have been made and tested. This technique seems to suit very 
well the operation of splicing optical waveguides without using a 
microscope and without fusing their tips together. The average results 
obtained in this way—power coupling efficiencies ranging from 80 per- 
cent to almost 90 percent, i.e., insertion losses under 1 dB and down to 
0.5 dB-compete successfully with the best results given by previous 
techniques. 

The facts, that rather mediocre end surfaces of the fibers can be 
accepted and that cladding diameters can differ at least by a few per- 
cent without affecting the performances, induce confidence in the 
future use of the sandwich technique for splicing multifiber cables in 
the field. It seems that all the operations could be extended to fiber 
tapes, dealing with each of them as a whole. A fairly precise alignment 
of the fibers would be required as the result of the tape manufacturing 
process, but the final alignment should be provided by a set of precisely 
embossed grooves on a plate. The fibers ought to be forced into the 
grooves by means of carefully designed mechanical tools, but without 
microscopic observation and without dealing with them on an indi- 
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vidual scale. The operations to be performed will depend on the ratio 
between the diameters of the fibers and their spacing in the tape, and 
on the precision of the preliminary alignment obtained in the manu- 
facturing process. The absence of stringent requirements on the quality 
of the end faces, proved by the available results, suggests that a simul- 
taneous and coplanar cut of all the fibers (for instance, scoring them by 
means of a razor blade while pulling along their axis) should be 
adequate. 
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B. 5. T. J. BRIEF 


Rain-Induced Differential Attenuation and Differential 
Phase Shift at Microwave Frequencies 


By J. A. MORRISON, M. J. CROSS, and T. 8. CHU 
(Manuscript received January 15, 1973) 


We give calculated results of differential attenuation and differential 
phase shift caused by rain, based on scattering of a plane electromag- 
netic wave by horizontally disposed oblate spheroidal raindrops. Two 
polarizations of the incident wave are considered, as depicted in Fig. 1. 
The factor e~**t has been suppressed. Subscripts I and II designate 
electric fields parallel and perpendicular to the plane containing the 
axis of symmetry of the raindrop and the direction of propagation of 
the incident wave, respectively. The angle between the direction of 
propagation and the axis of symmetry is denoted by a. For terrestrial 
microwave relay systems we take a = 7/2, but for satellite systems 
other values of a are of interest. 

The incident wave induces a transmitted field in the interior of the 
raindrop, and a scattered field. In the far field the quantities of pri- 
mary interest are the complex forward scattering functions! S8;(0) and 
S11(0). For the two polarizations considered, the polarization of the 
far scattered field is the same as that of the incident wave. However, 
when oblate raindrops are canted with respect to the vertical, the ver- 
tical and horizontal polarizations commonly used in radio communica- 
tions systems will be neither of the two aforementioned polarizations. 
Then cross-polarization coupling will occur as a result of the difference 
between S;(0) and S11(0). 

It is assumed that the ratio of minor to major axis of the oblate 
spheroidal raindrop depends linearly on the radius d@ (in cm) of the 
equivolumic spherical drop; specifically a/b = (1 — a). This relation- 
ship is similar to that used by Oguchi.2 The Laws and Parsons drop- 
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Fig. 1—Two polarizations of the incident wave. 


size distribution and Best’s raindrop terminal velocity are used.* Thus, 
for rain rates up to 150 mm/hour, there are 14 different drop sizes, 
a& = 0.025 (0.025) 0.35, to be considered. At 20°C, the refractive indices 
N = 5.581 + 2.8487 at 30 GHz and N = 6.859 + 2.7167 at 18.1 GHz 
were obtained from an elaborate fitting equation in a recently published 
survey‘ of available measured data. Since the calculations at 4 GHz 
were made at an earlier date, the value N = 8.77 + 0.9157 was taken 
from the older literature. 

Oguchi? has developed a perturbation theory for scattering from 
slightly eccentric spheroidal raindrops, but the first-order approxima- 
tion is expected to be inaccurate for the larger raindrops, which have 
larger eccentricity. Consequently, we have used a matching theory to 
obtain approximate nonperturbative solutions to the problem. Al- 
though the calculations reported on here are for oblate spheroidal rain- 
drops, the procedure may be used for axisymmetric raindrops which 
are not too nonspherical. Full details of the analytical and numerical 
procedures, together with tables of the computed values of S (0) and 
S311(0) for each drop size, will be given at a later date, but we outline 
the approach here. 

Spherical coordinates (r, 6, ¢) are chosen with polar axis along 
the axis of symmetry of the raindrop, and origin at the center, as 
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in Fig. 1. The scattered electric field is expanded in the form 
DZ (AmnM mn + OmaNmn), in terms of solutions of the vector wave equa- 
tion® satisfying the radiation condition. An analogous expansion is 
assumed for the transmitted field, in terms of vector wave functions 
which are finite at the origin. The complex coefficients in the expansions 
are determined approximately by satisfying the boundary conditions, 
namely the continuity of the tangential components of the total electric 
and magnetic fields across the surface of the raindrop, by a fitting pro- 
cedure. The index m denotes the order of the harmonic in the azimuthal 
angle y. The incident field may be expressed as a Fourier series in ¢, 
with coefficients depending on r and 6. Because of this, and the axial 
symmetry of the raindrop, the problem can be decomposed and the 
boundary fitting carried out independently for each m. 

A large number of points is chosen on the curve which is the inter- 
section of the boundary of the raindrop with the half-plane » = 0, 
0 < 6 S 7, and the coefficients are determined by requiring the bound- 
ary conditions to be satisfied at these points in the least squares sense. 
The advantage of using least squares fitting rather than collocation 
(in which the number of fitting points is equal to the number of un- 
known coefficients, which are then determined by solving a system of 
simultaneous linear equations) is that the boundary conditions are 
satisfied more accurately. At least twice as many fitting points as un- 
known coefficients were used in our calculations. The number of terms 
required to adequately satisfy the boundary conditions depends on 
both the frequency and drop size. For example, at 30 GHz it was neces- 
sary to take max m = 8 and max n = 238. In order to ensure the ac- 
curacy of the leading terms, it is necessary to take more terms in n 
than are really needed in the calculation of the far-field quantities such 
as S1(0) and S1;(0). A convergence test was carried out for each drop 
size, by increasing the upper limit of n in the sums by 2 and by 4. In 
a number of cases more than half the capacity of a Honeywell 6070 
computer was used. 

The rain-induced attenuation and phase shift are obtained from the 
forward scattering functions as follows :1 


2 
Aru = 0.434 — 2 Re Sr,11 (0) 9t (@) 43 /km 
Tv 


Py 11 


2 
—36 — » Im 81,11 (0) 0 (a) 92/km 
47? 
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Fig. 2—Rain-induced differential attenuation. 


where \ is the wavelength in centimeters, 9¢(@) is the number of rain- 
drops with mean radius @ per cubic meter, and the summation is taken 
over all drop sizes. For a = 90 degrees, the attenuation and phase shift 
have been calculated for various rain rates at 4, 18.1, and 30 GHz for 
polarizations I and II. The cases a = 50 degrees and 70 degrees have 
been calculated at 30 GHz. The differential attenuation Arr — Ar and 
the differential phase shift 3; — @; are summarized in Figs. 2 and 38. 
Full numerical results will be presented later. 

At 4 GHz the differential attenuation is negligibly small (0.036 
dB/km at 100 mm/hr), and therefore is not plotted in Fig. 2. How- 
ever, there is a differential phase shift of 5 deg/km at 100 mm/hr which 
can induce cross polarization as large as —10 dB over a long path of 
heavy rain; indeed, such a value was observed experimentally on one 
occasion by W. T. Barnett.® 

The calculated differential attenuations for a = 90 degrees in Fig. 2 
are considerably larger than those obtained by first-order perturbation 
theory. Taking into account the canting angle distribution,’ these 
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18.1 GHz, a =90°—— 


30 GHz 


DIFFERENTIAL PHASE SHIFT IN DEGREES/ km 





RAIN RATE IN mm /HR 
Fig. 3—Rain-induced differential phase shift. 


larger values of differential attenuation are needed to explain the mea- 
sured difference in attenuations for horizontal and vertical polariza- 
tions. ® 8 

At 18.1 GHz the differential attenuation and the differential phase 
shift will be comparable contributors to rain-induced depolarization. 
However, the differential phase shift remains small for heavy rain rates 
at 30 GHz, and hence the differential attenuation will be the dominat- 
ing cause of rain-induced depolarization at this frequency and pre- 
sumably also at higher frequencies. When a decreases from 90 degrees, 
both the differential attenuation and the differential phase shift be- 
come smaller as intuitively expected. 

The authors are indebted to D. C. Hogg for bringing this problem to 
their attention, to J. McKenna and N. L. Schryer for suggesting the 
matching and least squares fitting approaches, and for several helpful 
discussions in relation to these, to P. A. Businger whose least squares 
fitting subroutine was incorporated into the main program, to Mary 
Ann Gatto who took over the burdensome task of running the main 
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program, to Susan Hoffberg who wrote the program for calculating the 
first-order approximation for small eccentricities, and to Diane Vitello 
who performed the summation over the drop size distribution. The 
authors are particularly indebted to J. McKenna and D. C. Hogg for 
their continued encouragement throughout the lengthy course of this 
work. 


Note added in proof: 


A very recent paper by Oguchi,® in which similar calculations are 
carried out for a = 7/2 at 19.3 and 34.8 GHz, has come to our at- 
tention. He used collocation for the expansions in terms of spherical 
vector wave functions, and at 34.8 GHz he also used an expansion in 
terms of spheroidal wave functions. 
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